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ABSTRACT 



We examine a representative sample of 35 Seyfert 2 nuclei. Previous work has 
shown that nearly half (15) of these nuclei show the direct (but difficult-to-detect) 
spectroscopic signature at optical/near-UV wavelengths of the hot massive stars 
that power circum-nuclear starbursts. In the present paper we examine a va- 
riety of more-easily-measured quantities for this sample, such as the equivalent 
widths of strong absorption features, continuum colors, emission-line equivalent 
widths, emission line ratios and profiles, far-IR luminosities, and near-UV surface 
brightness. We compare the composite starburst + Seyfert 2 nuclei to "pure" 
Seyfert 2 nuclei, Starburst galaxies and normal galactic nuclei. Our goals are to 
verify whether the easily-measured properties of the composite nuclei are con- 
sistent with the expected impact of a starburst, and to investigate alternative 
less-demanding methods to infer the presence of starbursts in Seyfert 2 nuclei, 
applicable to larger or more distant samples. We show that starbursts do indeed 
leave clear and easily quantifiable imprints on the near-UV to optical continuum 
and emission line properties of Seyfert 2's. Composite starburst + Seyfert 2 sys- 
tems can be recognized by: (1) a strong "Featureless Continuum" (FC), which 
dilutes the Call K line from old stars in the host's bulge to an equivalent width 
W K < 10 A; (2) emission lines whose equivalent widths are intermediate between 
Starburst galaxies and "pure" Seyfert 2's; (3) relatively low excitation line-ratios, 
which indicate that part of the gas ionization in these Seyfert 2's (typically ~ 50% 
of H/3) is due to photoionization by OB stars; (4) large far IR luminosities (> 10 10 
L ); (5) High near-UV surface brightness (~ 10 3 L Q pc~ 2 ). These characteristics 
are all consistent with the expected impact of circum-nuclear starbursts on the 
observed properties of Seyfert 2's. Furtheremore, they offer alternative empir- 
ical diagnostics of the presence of circum-nuclear starbursts from a few easily 
measured quantities. 

Subject headings: galaxies:active - galaxies: Seyfert - galaxies: nuclei 



1. Introduction 

The role of starbursts in active galactic nuclei (AGN) is an issue which permeates the 
history of AGN literature. Theoretical scenarios for a "starburst- AGN connection" vary from 
those in which starbursts are a key piece of the AGN machinery, either through symbiotic 
or evolutionary processes (e.g., Perry & Dyson 1985; Sanders etal. 1988; Norman & Scoville 
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1988; Rees 1989; Daly 1990; Terlevich etal. 1992; Williams, Baker & Perry 1999), to those in 
which these are fundamentally decoupled phenomena, but which are likely to coexist merely 
because both are triggered by and live on gas fueling (e.g., Byrd etal. 1986; Byrd, Sundeluis 
& Valtonen 1987; Lin, Pringle & Rees 1988; Heller & Shlosman 1994; Hernquist & Mihos 
1995; Mihos & Hernquist 1996). Given the enormous breadth of possibilities in between 
these extremes, it is likely that progress in understanding this connection will be driven 
more by input from observations than from theoretical considerations. 

Observational work during the past decade has indeed substantially advanced our ability 
to recognize starbursts in the complex inner environment of AGN. Optical and UV studies 
of Seyfert 2 galaxies have been particularly enlightening in this respect. In a series of papers 
(Heckman etal. 1997; Cid Fernandes, Storchi-Bergmann & Schmitt 1998; Gonzalez Delgado 
etal. 1998; Storchi-Bergmann, Cid Fernandes & Schmitt 1998; Schmitt, Storchi-Bergmann 
& Cid Fernandes 1999; Storchi-Bergmann etal. 2000; Gonzalez Delgado, Heckman & Lei- 
therer 2001) we have detected unambiguous signatures of young massive stars within ~ 300 
pc of the nucleus in 30 to 50% of Seyfert 2's by means of high quality optical and, whenever 
possible, UV spectroscopy. This result fits well with the near IR studies of Oliva et al. (1995, 
1999), which find a comparable incidence of starbursts in Seyfert 2's through measurements 
of the stellar mass-to-light ratio. The starbursts in these composite starburst+AGN systems 
can make a significant contribution to the total luminosity output. In fact, studies specifi- 
cally aimed at candidate composite galaxies, as selected by emission line ratios intermediate 
between HII regions and AGN on classical diagnostic diagrams, demonstrate that their far 
IR and radio properties are dominated by star-forming activity at the 90% level (Hill et al. 
2001, 1999). Even the presence of compact radio cores, a more classical indicator of AGN 
activity (Condon etal. 1991; Sramek & Weedman 1986; Norris etal. 1990), can be accounted 
for by star- format ion (Hill etal. 2001; Kewley etal. 2000; Smith etal. 1998a,b; Lonsdale, 
Smith & Lonsdale 1993). Other recent results are covered in the review by Veilleux (2001) 
and the volume edited by Aretxaga, Kunth & Mugica (2001). 

Evidence is also steadily accumulating that associates AGN with star-formation on 
galactic scales. The ubiquity of super-massive black holes in the nuclei of normal galaxies 
in the local universe (Ho 1999), the proportionality between bulge and black-hole masses 
(Magorrian etal. 1998; Gebhart etal. 2000; Ferrarese & Merrit 2000), the link between 
Muack-hoie/ Mb u ig e and the age of the last major star-formation episode in the spheroid (Mer- 
rifield, Forbes & Terlevich 2000), all imply that the creation of these black holes and the 
ensuing QSO activity were an integral part of the formation of ellipticals and galactic bulges. 
This fossil evidence indirectly traces a much more active past, in which both copious star- 
formation and nuclear activity coexisted. Since the technology to directly study this high 
redshift era in detail is not available, we must guide our study of the starburst-AGN con- 
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nection by data gathered on nearby objects in which traces of this connection are caught in 
fraganti. This means refining our knowledge of circum-nuclear starbursts in Seyfert galaxies. 
In the optical-UV range, as hinted above, this is best accomplished studying type 2 Seyferts. 
Their favorable geometry, in which the blinding glare of the nucleus is blocked away by a 
dusty torus, facilitates the detection of features from circum-nuclear starbursts, thus making 
them the best-suited local laboratories to study the starburst-AGN connection. 

Advancing our empirical understanding of circum-nuclear starbursts in Seyfert galax- 
ies can be described as a 3-stage process: (I) Identifying starbursts in these systems; (II) 
characterizing the properties of the starburst (age and mass, or star-formation rate) and the 
AGN (black-hole mass, accretion rate); and (III) investigating possible connections between 
these properties. Once we get to this third stage fundamental questions can start to be ad- 
dressed. For instance: Is star-formation inextricably associated with AGN activity? Are the 
starburst and AGN powers related, say, by a M star b urst oc Muack-hoie proportionality between 
their masses? Or, analogously, is the star-formation rate connected to the rate of accretion 
onto the nuclear super-massive black-hole? Does the central engine evolve in parallel to the 
starburst around it? Much work remains to be done before tackling these key questions. 

The work on Seyfert 2's reviewed above pertains mostly to stage I of this process, as it 
essentially demonstrates how certain spectral features of massive stars can be used as sign- 
posts of starburst activity. Estimates of the recent history of star-formation in the circum- 
nuclear regions were also presented, thus touching the domain of stage II. These estimates 
come associated with known difficulties in accurately retrieving the detailed star-formation 
history of stellar systems (e.g., Leitherer 1999; Cid Fernandes etal. 2001; Kennicutt 1998), 
but they represent the current state-of-the-art regarding the characterization of the basic 
starburst properties. 

This paper deals mainly with stage I, that is, the identification of starbursts in Seyfert 
2's, using the data sets of Cid Fernandes etal. (1998) and Gonzalez Delgado etal. (2001, 
hereafter GD01) as a training set. This is probably the largest sample whose optical spectra 
have been systematically and carefully screened for starburst features, so it represents a 
good starting point. The analysis is geared towards verifying whether composite Seyfert 
2/starburst systems, those in which stellar features produced by young populations have 
been directly detected, exhibit other properties which differentiate them from systems not 
showing signs of starburst activity ("pure" Seyfert 2's). Unlike the detailed analysis reported 
by the previous work on these data, we here limit ourselves to easily measurable quantities. 
Specifically, we will use equivalent widths of strong absorption bands and continuum colors 
in the near UV, emission line equivalent widths, fluxes and line ratios, far IR luminosities 
and near UV surface brightness. The purpose of this analysis is two-fold: 
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(1) To verify whether the optical continuum and emission line properties of composite 
systems are consistent with the expected impact of starbursts upon these observables. 

(2) To investigate alternative ways to infer the presence of circum-nuclear starbursts in 
Seyfert 2's that may provide more straightforward diagnostics applicable to larger 
and/or more distant samples. 

The first of these points serves both as a consistency check and to give a sense of the 
role played by starbursts in the phenomenology and energetics of AGN. The relevance of 
the second point is that increasing the statistics is crucial to improve upon the still sketchy 
overall census of circum-nuclear starbursts in Seyfert 2's, to extend these studies to higher 
redshifts and to provide the raw material for further stage II and III studies. 

In §2 we describe the data for the Seyfert 2's and comparison samples used in this study. 
In §3 we deal with the problem of estimating the "Featureless Continuum" (FC) strength. 
We apply a method based on population synthesis, but which can be parameterized in terms 
of a few easily measurable quantities. Essentially, we are able to estimate the FC contribution 
to the optical light with little more than the equivalent width of the Call K absorption line, a 
method which may be of wide applicability in Seyfert 2 studies. In §4 we present an analysis 
of the emission line equivalent widths, line ratios and profiles, and discuss how they are 
affected by the presence of circum-nuclear starbursts. Emission line-FC correlations are also 
presented and discussed. In §5 we investigate how the composite systems behave in terms of 
far IR luminosity and how they fare in comparison to normal and interacting galaxies. In §6 
we compare the near UV surface brightness of Seyfert 2's in our sample to that of normal and 
Starburst galaxies. In §7 we collect the results obtained in a set of empirical criteria which 
may be used to diagnose the presence of circum-nuclear starbursts in Seyfert 2's. We also 
discuss the meaning of the composite/ "pure" Seyfert 2 classification and present tentative 
evidence of evolutionary effects in our sample. Section 8 summarizes our conclusions. 



2. The Database 

In this study we merge the Seyfert 2's in the southern sample studied by Cid Fernandes 
etal. (1998) with those in the northern sample of GD01, totaling 35 galaxies. Details of 
the observations, sample selection and previous analysis of these data can be found in the 
original papers and Heckman et al. (1997), Gonzalez Delgado et al. (1998), Storchi-Bergmann 
etal. (1998), Schmitt etal. (1999) and Storchi-Bergmann etal. (2000, hereafter SB00). We 
concentrate our analysis on the nuclear spectra, extracted through 1.2" x 2.1" apertures for 
the northern sample and 2" x 2" for the southern one. At the distances of our galaxies this 
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covers a region 60-860 pc in radius, with a median of 300 pc (H = 75 kms~ 1 Mpc~ 1 is 
adopted throughout this paper). We will also focus on the near UV to optical region of the 
spectrum, between 3500 and 5100 A. 

These 35 galaxies will hereafter be referred to as the "Seyfert 2 sample" . They constitute 
a good data set to investigate the starburst-AGN connection because the above papers have 
scrutinized these sources in search for signs of a starburst component. Specifically, the 
following signatures of a starburst were examined: (1) the presence of far-UV stellar wind 
lines (NV A1240, SHVA1397 and CIVA1549) due to O stars; (2) high order Balmer absorption 
lines of HI and Hel in the near UV; (3) the WR bump underneath HeIIA4686. The far-UV 
features are the cleanest signatures of recent star formation, but, unfortunately, there are 
only a handful of Seyfert 2's for which HST far-UV spectroscopy is feasible. All four galaxies 
in the sample for which we acquired such data (Mrk 477, NGC 5135, NGC 7130 and IC 
3639) have their far-UV spectrum dominated by young stars. For the remaining galaxies, 
detection of a starburst component relies on either high order Balmer lines, which originate 
in the photosphere of O, B and A stars, or the WR bump, a tracer of very recent (3-6 
Myr) star formation. Although more subtle than the far-UV features, these are reliable 
and more easily accessible indicators of the presence of a starburst. This is confirmed by 
the fact that either of these two features is also detected in the galaxies for which we have 
identified stellar wind lines in the far-UV. The least certain of these diagnostics is the WR 
bump, because independent evidence of the presence of WR stars is hard to obtain (GD01; 
Shaerer 2001). Nevertheless, this is the simplest interpretation of this feature, given the 
similarity to the feature sometimes seen in star- forming galaxies ("WR galaxies"), the lack 
of convincing alternatives and the simultaneous detection of other starburst features in some 
of the sources. 

The information to be used in this paper comprises fluxes and equivalent widths (W) of 
the [OII]AA3726,3729, HeIIA4686, H/3 and [OIII] A5007 emission lines; the W's of CaIIKA3933 
(W K ), CNA4200 (W C n) and G-bandA4301 (W G ) absorption bands, the continuum fluxes 
at 3660, 4020 and 4510 A, as well as the 12 and far IR luminosities as measured by 
IRAS. These are listed in Tables 1 and 2, along with heliocentric velocities extracted from 
NED 3 . Table 2 also lists the Ha/H/3 ratio for the Seyfert 2 sample, mostly compiled from the 
literature, since the Kitt Peak spectra do not cover the Ha region. These measurements come 
from spectra obtained with apertures larger than those used in our observations. For this 
reason, and also because of the intrinsic uncertainties associated with reddening corrections, 
we will concentrate on reddening insensitive diagnostics inasmuch as possible. Galaxies 



3 The NASA/IPAC Extragalactic Database (NED) is operated by the Jet Propulsion Laboratory, Cali- 
fornia Institute of Technology, under contract with the National Aeronautics and Space Administration. 



-7- 



which had hidden Seyfert 1 nuclei revealed through spectropolarimetry are indicated by the 
corresponding reference in the last column of Table 2. 

All absorption line Ws and continuum fluxes used in this paper were measured according 
to the system explained in Cid Fernandes etal. (1998). Emission line and continuum fluxes 
were corrected for Galactic extinction following the extinction law of Cardelli, Clayton & 
Mathis (1989, with R v = 3.1), and A B values from Schlegel, Finkbeiner & Davis (1998) as 
listed in NED. 



2.1. Composite and "pure" Seyfert 2's 

Fifteen of our Seyfert 2's have unambiguous evidence for starbursts as revealed by at 
least one of the 3 diagnostics outlined above. These are ESO 362-G8, IC 3639, Mrk 1, 
Mrk 78, Mrk 273, Mrk 477, Mrk 463E, Mrk 533, Mrk 1066, Mrk 1073, Mrk 1210, NGC 
5135, NGC 5643, NGC 7130 and NGC 7582 (observed prior to the type transition reported 
in Aretxaga etal. 1999), which will henceforth be tagged composite systems. The circum- 
nuclear starbursts in these galaxies cover a wide range of properties, from very young bursts 
still in the WR phase (e.g., Mrk 477, Mrk 1) to more mature systems in a "post-starburst" 
phase with deep high order Balmer absorption lines (e.g., Mrk 78, ESO 362-G8), and systems 
which have a mixture of both young and intermediate age stars (e.g., NGC 5135, NGC 7130). 

The remaining 20 sources will be labeled "pure" Seyfert 2's. This denomination is used 
simply to indicate that no clear signs of starburst activity have been detected in our previous 
studies. The spectra among our "pure" Seyfert 2's vary from galaxies whose optical light 
is entirely dominated by an old, red stellar population (e.g., NGC 1358 and NGC 5929) to 
galaxies which show clear signs of an UV-excess (e.g., Mrk 3 and Mrk 348). It is nearly 
impossible to tell on the basis of optical spectra alone whether this FC is associated with 
scattered light, a starburst or a combination of both. This point was discussed at length by 
SB00, who argue that as many as 30% of Seyfert 2's can fall in this ambiguous category. 
GD01 also finds that several galaxies are equally well described in terms of a scattered 
power-law or a starburst component, depending basically on what is chosen to represent the 
underlying old stellar population. Far-UV spectroscopy offers the best hope to break this 
degeneracy, but this is not practical with currently available instrumentation, so we will have 
to live with this ambiguity for some time. Indeed, this ambiguity pervades all discussions in 
this paper. 

We emphasize that the composite and "pure" categories are not meant as new elements 
in the AGN taxonomy. It is clear that this classification reflects a detectability effect, since 



Table 1. Seyfert 2 sample: Continuum, absorption and emission line properties 
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*Confirmcd starburst/Seyfert 2 composite 

Note. — Absorption line, continuum and emission line properties for the Seyfert 2 Sample, all extracted through small apertures (1.2" X 2.l" or 2" X 2"). All 
equivalent widths arc given in A and line fluxes in ergs -1 cm -2 (after correction for Galactic extinction). Empty slots correspond to weak lines, which were not 
measured. The flux calibration for NGC 5643, NGC 6300 and NGC 6890 are very uncertain due to weather conditions. 
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Table 2. Seyfert 2 sample: Additional data 



Galaxy 


A B,Gal 


cz 


Ha/H/3 


logil2 M 


log L FIR 




Rcf. 


Rcf. 


(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


CGCG420-015 


0.37 


8811 


5.5 


10.43 


9.76 


40.37 


dG92 




ESO 417-G6 


0.09 


4884 


4.7 


9.29 


8.45 


39.16 


S98 




BSO 362-G8* 


0.14 


4785 




9.34 


9.16 








Fairall 316 


0.60 


4950 


6.4 






38.80 


S98 




IC 1816 


0.12 


5080 


6.0 


9.75 


10.03 


39.98 


dG92 




IC 3639* 


0.30 


3275 


6.1 


9.98 


10.35 


39.79 


dG92 


H97 


IRAS 11215-2806 


0.31 


4047 


3.8 


9.56 


9.08 


38.88 


dG92 




MCG -05-27-013 


0.28 


7162 


5.0 


9.75 


9.99 


39.88 


dG92 




Mrk 1* 


0.26 


4780 


5.9 


10.47 


10.18 


40.03 


M94 




Mrk 3 


0.81 


4050 


6.6 


10.21 


10.18 


40.86 


M94 


MG90, T95 


Mrk 34 


0.04 


15140 


10.5 


10.33 


10.66 


40.69 


M94 




Mrk 78* 


0.15 


11137 


6.5 


10.34 


10.39 


40.24 


M94 




Mrk 273* 


0.04 


11326 


9.3 


10.62 


11.84 


39.88 


M94 




Mrk 348 


0.29 


4507 


6.0 


9.94 


9.84 


39.86 


M94 


MG90, T95 


Mrk 463E* 


0.13 


14895 


5.6 


11.19 


11.07 


41.25 


M94 


MG90, T95 


Mrk 477* 


0.05 


11332 


5.4 


10.35 


10.67 


41.09 


M94 


T95 


Mrk 533* 


0.25 


8713 


5.0 


10.85 


11.08 


40.54 


M94 


MG90, T95, H97 


Mrk 573 


0.10 


5174 


4.2 


9.66 


9.88 


40.16 


M94 


K94 


Mrk 607 


0.20 


2716 


4.9 


9.53 


9.63 


39.10 


D88 




Mrk 1066* 


0.57 


3605 


8.5 


9.90 


10.56 


40.06 


M94 




Mrk 1073* 


0.69 


6998 


6.3 


10.48 


11.04 


40.35 


M94 




Mrk 1210* 


0.13 


4046 


5.2 


10.05 


9.85 


40.28 


B99 


T95 


NGC 1068 


0.14 


1137 


4.5 


10.85 


10.79 


40.18 


M94 


AM85 


NGC 1358 


0.27 


4028 


3.4 


8.95 


9.34 


38.69 


M94 




NGC 1386 


0.05 


868 


4.9 


8.71 


9.09 


38.17 


M94 




NGC 2110 


1.62 


2335 


8.1 


9.42 


9.79 


39.28 


M94 




NGC 3081 


0.24 


2385 


4.5 






38.75 


M94 


MOO 


NGC 5135* 


0.26 


4112 


6.1 


10.17 


10.93 


39.99 


M94 




NGC 5643* 


0.73 


1199 


6.2 


9.34 


9.95 


38.77 


M94 




NGC 5929 


0.10 


2492 


6.2 


9.57 


10.21 


39.21 


M94 




NGC 6300 


0.42 


1110 




9.20 


9.80 








NGC 6890 


0.17 


2419 


4.2 


9.44 


9.87 


38.43 


M94 




NGC 7130* 


0.12 


4842 


6.6 


10.28 


11.05 


39.73 


S98 




NGC 7212 


0.31 


7984 


5.0 


10.24 


10.74 


40.61 


M94 


T95 


NGC 7582* 


0.06 


1575 


8.3 


9.74 


10.54 


39.11 


M94 





* Confirmed starburst/Scyfcrt 2 composite 

Note. — (1) Object name; (2) Galactic B band extinction (from NED); (3) cz in kms" 1 ; (4) Balmer decrement 
(compiled from the literature); (5) vL u at 12/j,m from IRAS, in Lq; (6) far IR. luminosity in Lq computed with the 
IRAS 60 and 100 fim fluxes using the formula in Sanders &i Mirabel (1996). When only upper limits were available 
for the IRAS fluxes we use half the quoted limit to compute fluxes. (7) H/3 luminosity in ergs -1 , corrected only for 
Galactic extinction; (8) reference for the Balmer decrement; (9) references regarding spcctropolarimctry detection of 
hidden Seyfert 1 nucleus. Code for references in columns (8) and (9): AM85- Antonucci & Miller (1985); B99- Bassani 
ctal. (1999); D88- Dahari &z Dc Robertis (1988); dG92 dc Grijp ct al. (1992); H97- Heisler, Lumsden & Bailey (1997); 
K94- Kay (1994); MOO- Moran ctal. (2000); M94- Mulchacy ctal. (1994); MG90- Miller & Goodrich (1990); S98- Schmitt 
(1998, PhD Thesis); T95- Tran (1995a, b). 
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we define composites as systems in which stars with ages of 10 6-8 yr have been conclusively 
detected. This definition therefore introduces an inevitable bias towards starbursts that are 
powerful in contrast with the other components of Seyfert 2 spectra, a bias which will become 
evident in several of the results presented in this paper. While it is likely that some "pure" 
sources eventually be "up-graded" to composites with deeper or different observations (e.g., 
when far-UV spectra are obtained, or when near-IR mass to light ratios are measured), it is 
unlikely that any of the composites will ever be reclassified as a "pure" Seyfert 2. Until such 
further notice, all galaxies not in the list of 15 composites above will be designated "pure" 
Seyfert 2's. 

The composite/ "pure" classification will be used throughout this paper as a guide to 
interpret our results, even though we shall not make explicit use of any of the properties on 
which this classification is based. Instead, we seek to find where composite sources fit into 
more general and easily accessible properties of Seyfert 2's. 

2.2. Comparison samples 

As comparison samples we draw data from three different observational studies. The first 
is the spectroscopy of 61 Starburst and Narrow Line Active Galaxies (NLAGN) by Storchi- 
Bergmann, Kinney & Challis (1995) and McQuade, Kinney & Calzetti (1995, hereafter 
collectively referred to as SBMQ95). The Starburst galaxies in their study will be referred to 
as the "Starburst sample" , which will be extensively used as a benchmark when gauging the 
results for our Seyfert 2's. UGCA410 was removed from our analysis because of apparently 
incompatible values of Wk and Wq, leaving a total of 36 Starburst galaxies. The remaining 24 
sources in SBMQ95 contain a mixture of Seyfert 2's, LINERs and starburst/AGN composites, 
so we give them a generic denomination of "NLAGN sample". In principle these galaxies 
serve as a test-bed to evaluate the efficacy of the diagnostics of compositeness developed 
from our Seyfert 2 sample. In practice, we will make limited use of this sample because 
of the mixed nature of its sources, which complicates the interpretation. Galaxies in the 
Starburst and NLAGN samples were observed with very large 10" x 20" apertures to match 
the IUE observations, sampling physical scales of typically ~1 kpc in radius at the distances 
of the sources. 

A third comparison sample used here is the "Merger sample" from the study of Liu 
& Kennicutt (1995a,b, hereafter LK95). This sample comprises 28 nuclear optical spectra 
integrated over a region corresponding to ~ 0.5-2 kpc in radius. The emission line properties 
of these galaxies are quite heterogeneous, comprising 5 Seyfert 2's, 9 HII nuclei and 14 
LINERs. Lutz, Veilleux & Genzel (1999) have used ISO data to show that IR-luminous 
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Fig. 1. — Aperture radii against the equivalent width of Call K for the Seyfert 2 (a), 
Starburst and NLAGN (b) and Merger (c) samples. Empty and filled circles in panel (a) 
indicate "pure" and composite Seyfert 2's respectively. Empty triangles in (b) correspond 
to NLAGN sources, whereas filled triangles correspond to Starburst galaxies. 
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galaxies with HII or LINER classifications are predominantly powered by starbursts, while 
dusty AGN are energetically significant in those classified as Seyfert 2's. Thus, we will use 
the Merger sample to empirically define the properties of powerful merger-induced dusty 
starbursts. 

The data which will be used from these samples consists of fluxes and W's of [Oil], H/3 
and [OIII]; Wk, Won and Wq, plus the continuum fluxes at 3660, 4020, and 4510 A and 
IRAS fluxes. Both the SBMQ95 and LK95 data offer more information than this, but we 
shall limit our analysis to the properties above for consistency with the Seyfert 2 sample. 

Fig. 1 shows the projected aperture radii, defined as r ap = (slit area/V) 1 / 2 , against 
Wk for the Seyfert 2, NLAGN, Starburst and Merger samples. Though the physical region 
sampled in the comparison samples extend well beyond that covered by the nuclear spectra 
of the Seyfert 2 sample, there is some overlap at small r ap which guarantees the fairness of 
the comparisons presented in this paper, and in any case aperture effects will be discussed 
whenever relevant. 

Fig. la also shows that the composite classification is not related to the physical size 
of the region sampled by our spectra. In principle one would expect that systems observed 
through larger apertures are more prone to exihibit starburst features due to the inclusion of 
light from off-nuclear star-forming regions. This effect has been quantified by Heckman et al. 
(1995) in their analysis of the aperture dependent properties of NGC 1068, whose 1 kpc ring 
of HII regions dominates the light for observations taken through kpc-scale apertures. Such 
rings are in fact rather common (e.g., Pogge 1989; Wilson etal. 1991; Storchi-Bergmann, 
Wilson & Baldwin 1996; Colina etal. 1997), and thus can play a role in defining a com- 
posite nature for more distant objects, as also discussed by Heckman etal (1995). If such 
aperture effects were pronounced in the Seyfert 2 sample we would expect composites to be 
concentrated towards large r ap , whereas Fig. la shows that composites and "pure" Seyfert 
2's are similarly distributed as a function of r ap . Furthermore, most of the apertures for 
the Seyfert 2 sample are much smaller than the 1 kpc typical of star-forming rings. This 
effect does not affect the Seyfert 2 sample, as our circum-nuclear starbursts are relatively 
compact (~ few hundred pc). These dimensions have in fact been confirmed by HST UV 
imaging for Mrk 477 (Heckman etal. 1997), NGC 5135, NGC 7130 and IC 3639 (Gonzalez 
Delgado etal. 1998), and are comparable to the inner star- forming structures identified in 
the Circinus galaxy (Wilson etal. 2000). A further example of how compact circum-nuclear 
starbursts around AGN can be is given in Shields etal. (2000). 

Another dataset that will be used in this paper is that of Jansen etal. (2000). This 
sample contains multicolor surface photometry for 198 galaxies, most of which classified as 
normal. We used their data to calculate the 3660 A surface brightness of these galaxies inside 



-13- 



radii of 250 and 500 pc, similar to the apertures used for our galaxies. These values are used 
to compare the UV surface brightness of Seyfert and normal galaxies. The 4 Seyfert 1 and 
1 BL Lac in their sample were excluded from our analysis. 

3. Analysis of the Starlight and the Featureless Continuum in Seyfert 2's 

Determining the FC contribution to the optical spectrum of Seyfert 2's is an important 
but difficult task, for which there is no definitive recipe. The traditional and still most widely 
employed method, which dates back to at least Koski (1978), is to remove the contribution 
of the old red stars from the host's bulge by adopting an elliptical or early type spiral normal 
galaxy as a spectral template. The spectral decomposition is accomplished by assuming a 
F v oc v~ a power-law for the "non-stellar" FC, whose strength is then adjusted to minimize 
the residuals, yielding typically a = 1-2. Naturally, this procedure precludes the detection 
of young stellar populations, whose most conspicuous optical absorption features are masked 
by emission lines even in Starburst galaxies. Furthermore, there is no simple way to tell 
apart a power law FC from that produced by a reddened starburst (e.g., Cid Fernandes & 
Terlevich 1995). Only upon careful scrutiny of high quality spectra is it possible to discern 
the presence of a starburst from its high order Balmer lines or, sometimes, the WR bump 
(see discussion in GD01). Therefore, the FC strength will in general include a contribution 
from young stars as well as a bona fide non-stellar FC, composed of scattered light from the 
hidden nucleus plus nebular continuum. 

"Featureless" is of course a misnomer for a spectral component at least partly associated 
with starlight. We nonetheless retain this terminology both for historical reasons and because 
to first order a young stellar population does produce weak spectral features in the optical 
range. When referring to scattered light from the hidden nucleus we will use the term 'FC1' 
introduced by Tran (1995c). As argued by Cid Fernandes & Terlevich (1995), Heckman 
etal. (1995) and Tran (1995c), strong FC Seyfert 2's must also contain another source of 
FC photons, FC2' in Tran's notation, to account for the absence of conspicuous broad 
lines and the larger polarizations observed in the reflected broad lines than in the reflected 
nuclear continuum. At least for the composite systems, there is little doubt that FC2 is 
associated with circum- nuclear starbursts (see GD01). Nebular continuum and emission 
from the scattering region itself (Tran 1995c) are other possible contributors to FC2. 

In this paper we explore a different way to determine the FC strength. We decompose 
the total light into a base of 12 simple stellar population components represented by star 
clusters spanning the 10 6 to 10 10 yr age range and different metallicities (Schmidt et al. 1991; 
Bica 1988), plus a z/~ L5 power-law to represent FC1 and any other AGN component. This is 
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done by means of a population synthesis analysis similar to that employed by Schmitt et al. 
(1999), but performed with the probabilistic formalism described in Cid Fernandes et al. 
(2001), extended to include the power-law component. Note that this method is analogous, 
but not technically equivalent to a full spectral decomposition as only a few W's and colors 
are actually synthesized. As input data we use the W s of Call K, CN and G-band, plus the 
-F3660/-F4020 and F4510 /F4020 colors. The errors on these quantities were set to 0.5 A for Wk 
and W G , 1 A for Won and 0.05 for the colors; these errors are consistent with the quality 
of the spectra, and in any case they do not affect our conclusions. The output is a 13-D 
population vector x which contains the expected values of the fractional contribution of each 
component to the total light at a reference normalization wavelength. In what follows we 
shall use 4861 A as our reference wavelength unless noted otherwise. Uncertainties in x are 
also provided by the code. 

It is clear from the outset that with so little input information we will not be able to 
recover all 13 components of x accurately, but, as discussed by Cid Fernandes etal. (2001), 
we may rightfully hope to achieve a coarse yet useful description of the population mixture 
by grouping similar elements of the population vector. In what follows we shall make use of 
only three grouped x components: Xold, made up from the sum of all base components with 
age > 1 Gyr, xint, corresponding to the "post-starburst" 10 s yr intermediate age bin, and, 
most importantly for our purposes, Xfc, containing the total contribution of < 10 7 yr stars 
plus the power-law component. These quantities are linked by the xold + xint + xfc = 1 
normalization constraint, so in practice our method characterizes the spectral mixture in 
terms of 2 parameters. This bi-parametric description of the data is in many ways analogous 
to a Principal Component Analysis (Sodre & Cuevas 1997; Rodrigues-Lacerda 2001), except 
that by construction we have a priori knowledge of the physical meaning of our components. 

Note that we do not attempt to disentangle the starburst and power-law components 
explicitly; these are merged in xpc — x ys + %pl-, where xys an d xpl stand for the light 
fractions due to young (< 10 7 yr) stars and the power-law respectively. As shown below, 
both the FC strength and the emission line properties give us strong hints as to which of 
these two components dominates xpci thus helping to break this degeneracy. 

3.1. Results of the Synthesis 

The results of this analysis are listed in Table 3 in order of decreasing xfc- A first 
noticeable result in this Table is that composites are concentrated towards the top of the 
list: All 10 xpc > 29% sources are confirmed Seyfert 2/starburst composites, and 12 of the 
15 composites are among the 16 sources with xpc > 15% at 4861 A. The identification of 
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Table 3. Synthesis Results 



Galaxy 


%OLD 


[%] 


XINT 


[%] 


XFC [%] 


(1) 


(2) 




(3) 




(4) 


TV /T 1 A T~7-k 

Mrk 477 


1U.0 


_L 
ZL 


3.4 


9 9. 
Z.O 


ZL 


2.0 


8o. o ± 4.2 


Mrk 400.CJ 


97 7 


4_ 

zl 


4.3 


O.O 


4_ 

ztz 


2.5 


OO.O ± 4.U 


mo o cioc* 
InCtL- 51o5 


97 9 
Z 1 .Z 


_L 
ZL 


7.3 


OR Q 

zo.y 


4_ 
ZL 


6.2 


4o.y ± o.o 


MO O 7ion* 
JNCxL- / loU 


OO.O 


_U 

ZL 


8.0 


ZO.O 


4_ 
ZL 


6.1 


A A 1 _l_ c n 

4U.1 ± o.y 


Mrk 5oo 


40. Q 


4- 

ZTZ 


( . 1 


9^ 3 
zo.o 


4- 

ZL 


6.1 


33.8 i 5.5 


Mrk lOoo 


01.0 


4- 

ZL 


6.8 


1 « 3 

lO.O 


4- 

ZL 


5.4 


32.3 ± 4.9 




OO.O 


_L 
ZL 


7.8 


OO.O 


4_ 

ZL 


6.5 


on a _l c T 
oU.o ± 5.7 


Mrk l 


Pil 1 
Oo. 1 


4_ 

ZL 


O.O 


7 1 


4- 

ZL 


O.O 


OO Q J- a n 
zy.o it 4.U 


Mrk 1U { 6 


OO.O 


4- 

ZL 


6.2 


1 ^ K 
10.0 


4- 

ZL 


5.1 


29.1 i 4.6 


l\ /I-wlr 1 O 1 n* 

Mrk izlU 


DO. 1 


_L 
ZL 


A *7 


^ Q 

o.y 


4- 

ZL 


3.1 


on r> I o f 

zy.u ± o.o 


iNtjL/ lUoo 


fi^ 9 
OO.Z 


4_ 
ZL 


0.0 


ii.i 


4- 

ZL 


4.3 


OK ^7 _L a n 

25.7 ± 4.U 


Mrk o 


77 S 
1 / .O 


4_ 
ZL 


A. ( 


i i 
i . i 


4- 

ZL 


0.8 


21.1 i 2.0 


MOO ^70 1 O 

IN^jrL/ ( AY A 


77 ^ 
I f .o 


4- 

ZL 


3.6 


9 n 

z.u 


_|_ 

ZL 


1.4 


on o _i_ o a. 
2U.o ± z.O 


TV /T 1 iT"71* 

Mrk 27 o 


*± i .O 


_L 
ZL 


7.1 


o^t.o 


4- 

ZL 


6.2 


17.9 ± 4.7 


1\ /Tv.1, OA 

Mrk o4 


78.1 


± 


4.6 


4.1 


± 


2.5 


17.7 i 3.1 


TO QftQO* 


69.9 


± 


O.O 


14.4 


± 


O.U 


10. / ±4.1 


Mrk o4o 


87.2 


± 


A A 

4.4 


3.2 


± 


2.0 


9.7 ± 2.6 


NGC 2110 


87.8 


± 


3.4 


2.5 


± 


1.6 


9.6 ± 2.0 


IC 1816 


84.5 


± 


4.8 


7.2 


± 


3.3 


8.3 ± 2.6 


NGC 5643* 


73.2 


± 


7.2 


19.5 


± 


5.3 


7.3 ± 3.2 


Mrk 78* 


72.0 


± 


6.5 


21.2 


± 


4.9 


6.8 ± 2.9 


NGC 3081 


90.1 


± 


4.3 


3.4 


± 


2.0 


6.4 ± 2.1 


Mrk 573 


90.9 


± 


4.3 


3.4 


± 


2.0 


5.6 ± 1.9 


NGC 5929 


91.7 


± 


3.7 


2.7 


± 


1.7 


5.6 ± 1.9 


NGC 6890 


93.1 


± 


4.5 


1.9 


± 


1.4 


5.0 ± 2.2 


NGC 1386 


89.9 


± 


5.4 


5.1 


± 


2.7 


5.0 ± 2.1 


Mrk 607 


90.4 


± 


4.4 


4.7 


± 


2.4 


4.9 ± 1.8 


ESO 362-G8* 


56.6 


± 


5.9 


38.5 


± 


4.9 


4.9 ± 2.4 


CGCG 420-015 


91.7 


± 


5.0 


4.1 


± 


2.3 


4.2 ± 1.8 


MCG -05-27-013 


93.9 


± 


3.7 


1.9 


± 


1.3 


4.2 ± 1.6 


IRAS 11215-2806 


90.0 


± 


5.7 


5.9 


± 


3.0 


4.0 ± 1.9 


ESO 417-G6 


92.7 


± 


4.6 


4.0 


± 


2.2 


3.3 ± 1.5 


Fairail 316 


96.6 


± 


3.8 


1.5 


± 


1.0 


1.9 ± 0.9 


NGC 6300 


96.5 


± 


4.6 


1.7 


± 


1.2 


1.8 ± 1.0 


NGC 1358 


97.6 


± 


3.8 


1.2 


± 


0.8 


1.2 ± 0.7 



*Confirmcd Starburst/Seyfert 2 composites 

Note. — (1) Object name; (2-4) fractions of the total light 
at 4861 A due to old stars (xold), 100 Myr stars (xint)\ 
and < 10 Myr old stars plus a power-law FC (xfc)- 

Galaxies are organized in order of descending FC strength. 
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Fig. 2. — Results of the synthesis analysis, condensed into a bi-parametric xpc an d £/att 
representation. Dotted lines trace lines of constant xold, as labeled, (a) Sources from the 
Seyfert 2 sample, with filled circles indicating the starburst /Seyfert 2 composites. The empty 
circle among the xold < 75% zone otherwise exclusively occupied by composites is NGC 
1068. As xold increases it becomes progressively difficult to detect starburst features in 
Seyfert 2's, explaining the clear separation between composites and "pure" Seyfert 2's in 
this diagram, (b) Results for the Starburst (filled triangles) , NLAGN (empty triangles) and 
Merger (filled squares) samples. 
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xfc with xys (i-e., xys 3> ^pl) is thus safe for these objects. 

An even more remarkable result is that all 15 composites but only one "pure" Seyfert 2 
(NGC 1068) have x ld < 75%. This is illustrated in Fig. 2a, where we condense the results 
of the synthesis in the xfc x %int plane, with dotted lines indicating lines of constant 
Xold = 1 — Xint — xfc- The diagram segregates composite from "pure" systems very 
effectively, and it does so in agreement with what one would expect: composites have larger 
proportions of intermediate age and young stars, the latter being included in xpc- 

As a check on our method, we have applied it to the spectra in the comparisons sam- 
ples, which, for consistency, were processed in exactly the same way as our Seyfert 2's. If 
our method works, we would expect bona fide Starburst galaxies to be located within the 
Xold < 75% zone, probably more towards even lower values of Xold, overlapping with our 
more extreme composites. This is confirmed in Fig. 2b, which shows the synthesis results for 
the Starburst, Merger and NLAGN samples. There is not a single Starburst galaxy (filled tri- 
angles) for which xold > 65%, whereas most of the NLAGN are located there. Furthermore, 
most of the NLAGN which intrude into the Starburst region are known starburst /Seyfert 2 
composites! The four empty triangles in the x ld < 50% zone, for instance, are Mrk 477, 
NGC 5135, NGC 7130 (all also in our Seyfert 2 sample) and NGC 7496, whose composite 
nature is discussed in Veron, Gongalves & Veron-Cetty (1997). The 50% < xold < 75% 
zone also contains known composites, both from our sample (IC 3639 and NGC 7582), and 
from independent work, such as NGC 4569 (Keel 1996; Maoz etal. 1998), NGC 1672 (Veron 
etal. 1997) and NGC 6221 (Levenson etal. 2001b). 

As expected, most of the sources in the Merger sample (squares in Fig. 2b) behave 
like starbursts in our stellar population diagram. Note that our synthesis indicates a strong 
contribution of ~ 10 8 yr populations for most sources in this sample. This is consistent with 
the results of LK95, who find a high incidence of systems in a post-starburst phase based on 
the detection of the Balmer absorption series. Those mergers which are dominated by old 
stars (xold > 50%) are predominantly LINERs (e.g., NGC 942, NGC 3656 and 3C293). 

It is thus clear that, despite its limitations, our bi-parametric description of the data 
provides a very efficient empirical diagnostic of compositeness in Seyfert 2's. Since our 
operational definition of composites is based on the detectability of starburst features, which 
is certainly facilitated when they make a strong contribution to the continuum, the result 
that sources with large xpc and/or xj^t are mostly composites is somewhat redundant. 
The importance of this result is twofold. First, it shows a striking empirical similarity 
between the bona fide Starburst galaxies (whose optical/near- UV continuum is definitely 
produced by young and intermediate-age stars) and the starburst /Seyfert 2 composites. 
Second, it demonstrates an excellent consistency between the two techniques, which are 
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based on different observables. The synthesis process did not use any of the information 
used in the original identification of starbursts in these composites, namely, far-UV stellar 
wind features, high order Balmer lines, and/or the WR bump, all of which are much more 
subtle or hard to obtain than the observables used as input for the synthesis code. 

3.2. Contrast and Evolution 

The separation between composites and "pure" Seyfert 2's in Fig. 2a strongly suggests 
that, as anticipated in §2.1, the composite/ "pure" classification is mostly driven by a contrast 
effect: Systems with circum- nuclear starbursts residing in galaxies where the old stellar 
population dominates the spectrum would simply not be recognized as composite systems. 
This contrast effect is nicely quantified by the value of xold, which can therefore be read 
as a measure of the difficulty to detect any non-trivial spectral component and to interpret 
its origin. Fig. 2 shows that our current techniques establish a threshold of x OLD ~ 75%, 
beyond which we are not able to recognize circum-nuclear starbursts. For instance, if we 
scale down the circum-nuclear starburst in IC 3639 by a factor > 2, it would move from 
xold = 70% to more than 82%, where, judging by Fig. 2, we would classify it as a "pure" 
Seyfert 2. It is therefore likely that at least some "pure" Seyfert 2's are just that: composites 
with weak starbursts. 

The presence of strong HI emission lines introduces another contrast effect, since they 
dilute the high order Balmer absorption lines, thus hindering the use of this diagnostic of 
starbursts. Therefore the combination of strong emission lines and a large xold is highly 
unfavorable to the detection of circum-nuclear- starbursts. We will return to this issue in 
§4.1.1, after we analyze the emission line data. 

Besides the contrast between the starburst and the underlying old stellar population, 
evolution is the main property defining the location of a source in Fig. 2. One can easily 
imagine an evolutionary sequence which starts at some large value of xpc in the bottom right 
part of the diagram, with massive young starts dominating the spectrum and no intermediate 
age population, and gradually moves towards large xint/xfc (i-e., to the top left of the 
diagram) on ~ 10 8 yr time scales. The exact shape of evolutionary tracks in our xfc x 
xint diagram is dictated by the detailed star-formation history, but for the purposes of the 
qualitative discussion below one can imagine that such tracks broadly follow the dotted lines 
traced in Fig. 2 up to a few hundred Myr, and then collapse to the origin, after the starburst 
ends and its stars move into our old population bin. 

Here too we find a remarkable agreement between the present synthesis analysis and the 
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results of our previous work on the composites in the Seyfert 2 sample. The four composites 
with larger xfc/ x int are Mrk 477, Mrk 463E, Mrk 1210 and Mrk 1 (going upwards in 
Fig. 2a). These are precisely the 4 systems where we have identified WR stars (Heckman 
etal. 1997; Storchi-Bergmann etal. 1998; GD01), signposts of very recent star-formation. 
They therefore rank as the youngest circum-nuclear starbursts among our composites, in 
agreement with their location in Fig. 2a. At the other extreme, the three composites with 
x FC < 10% are NGC 5643, Mrk 78 and ESO 362-G8, all of which present spectroscopic 
signatures of a dominant 'post-starburst' population (Schmitt etal. 1999; SB00; GD01), and 
rank as our oldest composites, in agreement with the large Xint but small Xfc we obtain. 
Given this agreement at the young and old ends of a starburst evolutionary sequence, we 
would expect that systems located at xfc/%int values in between these extremes present 
a mixture of young and intermediate age stars. This is indeed the case. The cluster of 4 
composites towards the top right in Fig. 2a, for instance, contains NGC 5135, NGC 7130, 
NGC 7582 and Mrk 533, all of which reveal spectroscopic signatures of this mixture (SB00; 
GD01), such as pronounced high order Balmer absorption lines simultaneous with far-UV 
stellar wind lines. 

What makes this agreement remarkable is the fact that while the characterization of 
the starburst component in our earlier work involved a detailed spectroscopic analysis, the 
synthesis performed here used just three Ws of metallic absorption bands plus a couple of 
near-UV colors. Our quicker and cheaper method is therefore not only able to recognize 
composite systems from a handful of easily measured observables, but it allows us to go one 
step further to provide a rough description of evolutionary state of the starburst component. 

We will return to the issue of evolution in §7. Before that, we have to deal with 
the contribution of the non-stellar FC in Xpc- This issue was deliberately omitted from 
the discussion above because we do not have a recipe to separate the stellar and non-stellar 
parts of the FC, neither do we know how or whether the non-stellar FC evolves on time-scales 
comparable to the starburst lifetime. These difficulties limit considerations about evolution 
to composites alone, for which we know that the starburst component is the dominant 
contributor to the FC. 



3.3. The FC component 

Our inability to disentangle the starburst and non-stellar components of the FC is 
currently the major obstacle to a full characterization of the starburst- AGN connection in 
Seyfert 2's. Unfortunately, these are also the most relevant spectral components in Seyfert 
2's. The starburst portion of x FC traces the ongoing star-formation, being therefore asso- 
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ciated with young massive stars which may have a significant impact on the ionization of 
the gas in the circum-nuclear environment. On the other hand, the non-stellar component 
(which is presumably mainly scattered light, i.e., FC1) is the only continuum feature associ- 
ated with the AGN in these galaxies. For these reasons we concentrate on the FC component 
throughout most of the rest of this paper. Whereas little can be said about the starburst 
and AGN shares of xfc in "pure" Seyfert 2's, composite systems have their FC dominated 
by the starburst (GD01; SB00), and hence provide a useful guide as to what can be deduced 
from the FC strength alone. 

3.3.1. Strong FC sources are Seyfert 2/ starburst composites 

Table 3 shows that strong FC sources (xfc ^ 30%) are likely to be composite systems. 
Within the context of the unified model, it is not surprising to find that a starburst com- 
ponent dominates over FC1 whenever the total FC is strong, since scattered light cannot 
exceed a fraction of ~ 30%, otherwise reflected broad lines should become easily discernible 
in the direct spectrum and the galaxy would no longer be classified as a type 2 Seyfert. For a 
typical W of ~ 100 A for broad H/3 in Seyfert l's (Goodrich 1989; Binette, Fosbury & Parker 
1993), more than 30% scattered light would imply a broad H/3 stronger than 30 A in the 
direct total spectrum. Cid Fernandes & Terlevich (1995) argue that broad lines should be 
seen even for smaller xfci, but we figure 30% is a reasonable limit in practice because NGC 
1068 has xfc = 26%, and broad lines have not been conclusively detected there without 
the hindsight of spectropolarimetry to guide the eye (Miller & Goodrich 1990; Tran 1995b; 
Malkan & Filippenko 1983). 

In support of this conclusion, we note that several cases of strong FC sources reported 
in the literature turned out to harbor starbursts not accounted for in the classical starlight 
template decomposition. For instance, Miller & Goodrich (1990) find xfc = 100% in Mrk 
463E, while GD01 suggest it actually contains a powerful starburst in the WR phase. The 
same can be said about Mrk 477, for which Tran (1995a,b) finds x F c — 59% at 5500 A, and 
Mrk 1066, for which Miller & Goodrich (1990) find a 72% FC at 5300 A. Even in Cygnus 
A, whose blue FC has remained a mystery for more than a decade (e.g., Goodrich & Miller 
1989; Tadhunter, Scarrott & Rolph 1990), there now seems to be evidence for young stars 
(Fosbury etal. 1999). Indeed, inspection of the Keck spectrum in Ogle etal. (1997) even 
indicates the presence of a possible WR bump! This is consistent with our conclusion, since 
estimates of the FC strength for Cygnus A place it above the 30% level (Osterbrock 1983). 

This high "success" -rate suggests that one can predict the existence of circum-nuclear 
starbursts based only on the FC strength. Despite the non uniformity of FC estimation 
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methods, one can be fairly confident that a search for sources with xfc > 30% reported in 
the literature should identify other composite systems. Kay (1994), for instance, reports the 
following Seyfert 2's with an FC stronger than 30% at 4400 A: Mrk 266SW (30%), Mrk 1388 
(45%), NGC 591 (46%), NGC 1410 (38%), NGC 1685 (38%), NGC 4922B (53%), NGC 7319 
(35%), NGC 7682 (30%). Koski (1978) finds a 34% FC at 4861 A for 3C184.1 and 35% for 
Mrk 198. Evidence for young and intermediate age stars in the spectra of Mrk266SW was 
reported by Wang etal. (1997). According to Gongalves, Veron-Cetty & Veron (1999), this 
galaxy also has a composite emission line spectrum, with emission line ratios intermediate 
between those of Seyferts and Starbursts. The same description is given for the interacting 
galaxy NGC 4922B by Alonso-Herrero et al. (1999). We could not find information regarding 
the compositeness or otherwise of the remaining galaxies. Our prediction is that most of 
them should present detectable signatures of circum-nuclear starburst activity. 



3.3.2. Moderate and weak FC sources are ambiguous 

While all xfc ^ 30% sources are composites, the converse is not true. For IC 3639, a 
galaxy whose composite nature has been conclusively established by both optical and UV 
spectroscopy, we find a 16% FC. The range between 15 and 30% also contains 4 "pure" 
Seyfert 2 nuclei: NGC 1068, Mrk 3, NGC 7212 and Mrk 34. 

NGC 1068 is sometimes listed as an example of a starburst-AGN link because of its 
bright ring of HII regions at ~ 1 kpc from the nucleus, which contributes roughly half of the 
bolometric (IR) luminosity of this prototype Seyfert 2 (Telesco et al. 1984; Lester et al. 1987). 
If observed from further away, the ring would dominate the optical spectrum, imprinting 
signatures which would cast distant NGC 1068's in the composite category (Heckman etal. 
1995; Colina etal. 1997), but, as argued in §2, this is not the effect behind the starburst 
features identified in our composites (Fig. 1). The spectrum of NGC 1068 analyzed here 
pertains to a much smaller region, corresponding to r ap = 66 pc. The only indications of star- 
formation within this central region are the strong Call triplet at ~ 8500 A (Terlevich, Diaz 
& Terlevich 1990) and the small M/L ratio at 1.6 fan (Oliva etal. 1995, 1999), both of which 
point to an intermediate- age population, qualitatively consistent with the xint = 11% found 
in our synthesis. There does not seem to be substantial ongoing star-formation associated 
with this "post-starburst" population, since neither far-UV (Caganof etal. 1991) nor optical 
(GD01; Miller & Antonucci 1983) spectroscopy reveal signs of young stars close to the 
nucleus. Furthermore, from the work of Antonucci & Miller (1985), Miller & Goodrich 
(1990) and Tran (1995c) we know that NGC 1068 does not suffer from the "FC2-syndrome" 
(lower polarization in the continuum than in the broad lines), so most, if not all, of its nuclear 
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FC is indeed FC1. Since our xfc estimate is in fair agreement with previous determinations 
(e.g., the 22 % at 4600 A found by Miller & Antonucci 1983 and the 16% at 5500 A derived 
by Tran 1995), there is little doubt that xpci 3> x ys m the central region of NGC 1068. We 
therefore keep NGC 1068 in the "pure" Seyfert 2 category. 

The nature of the FC is not so clear for the 3 other moderate FC "pure" Seyfert 2's 
in our sample: Mrk 3, NGC 7212 and Mrk 34. In their detailed analysis of the same data 
used here, GD01 find that the dilution of the metal absorption lines in these galaxies is 
better modeled by an elliptical galaxy plus power-law spectral decomposition than using an 
off-nuclear template, which, together with the lack of starburst features, lead them to favor a 
scattered light origin for their FC. The extended blue continuum aligned with the radio axis 
and ionization cone in Mrk 3 (Pogge & De Robertis 1993) supports this idea. Kotilainen 
& Ward (1997) find a similar feature in NGC 7212, though not aligned along the [OIII] 
emission. These indications of a FC1 component are confirmed by the spectropolarimetry 
observations of Tran (1995a,b,c), which revealed their hidden Seyfert 1 nuclei. However, 
Tran also finds that Mrk 3 and NGC 7212 suffer from acute "FC2-itus"! According to his 
analysis, only 4% (Mrk 3) and 2% (NGC 7212) of the flux at 5500 A is attributable to FC1; 
the rest of their FC is 'FC2' ! The total FC = FC1 + FC2 fractions at 5500 A found by 
Tran are 12% for Mrk 3 and 17% for NGC 7212, which are in very good agreement with our 
estimates of xfc when translated to the same wavelength. Therefore, from the point of view 
of their polarization spectra, most of their FC still have to be explained. 

Below xfc ~ 30% we therefore enter a "gray zone": the power-law/starburst degener- 
acy sets in and the identity of the FC becomes ambiguous. Extra information (spectropo- 
larimetry, UV spectroscopy, imaging) may help disentangling these two components in a few 
sources, but the situation becomes even fuzzier for weak FC sources (x FC < 15%). The 
case of Mrk 348 is illustrative in this respect. Using an elliptical galaxy template, Tran 
(1995a,b) finds a 27% FC at 5500 A, which, when combined with his spectropolarimetry 
data, propagates to the conclusion that 22% of the total light is associated with something 
else, ie., FC2. Storchi-Bergmann etal. (1998), on the other hand, find a much weaker FC by 
using an off-nuclear spectrum as a starlight template, a conclusion corroborated by GD01 
on the basis of independent observations. In fact, they find that the nuclear spectrum is well 
matched with no FC at all, but, guided by the results of Tran, they favor a model in which 
5% of the light comes from the FC1 component detected via spectropolarimetry. 

These inconsistencies undoubtedly emerge due to the intrinsic weakness of the FC, 
which, combined with the dispersion of stellar populations in active galaxies (Cid Fernandes 
etal. 1998; Schmitt etal. 1999; Boisson etal. 2000; GD01), boosts the differences between 
different FC estimation methods. Whereas in Mrk 348 the combination of spectropolarimetry 
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and long-slit spectra allowed considerations about the nature of its FC, we do not have such 
information for most sources below Mrk 348 in Table 3, and even if such data existed, it 
is clear that any attempt to perform a detailed spectral decomposition would be highly 
uncertain, since it is simply unrealistic to expect accuracies better than ~ ±5% in any 
xfc determination method. From the point of view of this paper, we consider all such 
sources ambiguous. We hope to shed light on the nature of these moderate and weak FC 
Seyfert 2's by considering information other than that contained in the optical continuum 
and absorption lines (§4). 



3.4. Empirical calibration of the FC strength 

As in any other study dealing with optical spectra of Seyfert 2's, the FC strength plays 
a major role throughout this paper. The results above already indicate that, regardless of 
the ambiguities involved in the interpretation of xpc, our method was able to recover a 
meaningful component of the near UV-optical continuum, a conclusion which will become 
stronger in the analysis that follows. This in itself is an important result, since we based 
our FC estimation on just a handful of easily measurable quantities. Though the synthesis 
process is rather elaborate, we may use its results to derive an a posteriori calibration of 
xpc in terms of the observables it is based upon. In fact, as shown in Fig. 3, we find that 
the expression 
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„„= -0.33 (^j +0.52 ^ + 0.89^-1.04^j^]-0.08. (1) 

where Wk is in A, recovers xfc within ±3% for all galaxies. Since the measurements of 
both Wk and the -F3 660 /F4020 color do not require large S/N spectra, this calibration can be 
applied to data sets not meeting the S/N > 25 standard of our sample (which was necessary 
to identify weak starburst features). This provides a much more straight-forward and well 
defined way to estimate xfc- 

Although old stars are of no direct interest in the starburst- AGN connection, we have 
seen that they play a central role in defining the detectability of circum-nuclear starbursts, 
so it is useful to have an equation analogous to (1) for x OLD . We find that 

xold = 1-92 - 0.98 H^f\ 2 + 0.05 (2) 
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Fig. 3. — (a) Relation between the FC strength at 4861 A and the equivalent width of Call 
K for the Seyfert 2 sample. Composites are plotted with filled circles, (b) Comparison of the 
FC strengths derived from the synthesis and from the empirical fit of xpc as a function of 
Wk and a near UV color (equation 1), which offers an easy way of estimating xpc- Starburst 
galaxies from the SBMQ95 data are also included in this plot (filled triangles), but not in 
the fit. 
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reproduces the synthetic x ld to better than ±5% for all galaxies. 

Equations (1) and (2) also do a good job for the Merger, NLAGN and Starburst samples. 
This is illustrated in Fig. 3b for the galaxies in the Starburst sample. The difference, of 
course, is that for these bona fide Starbursts the FC is undoubtedly produced by young 
stars. 

A corollary of x FC and x OLD being so closely related to Wk is that this quantity can 
be used as a first order measure of the stellar population mix and the spectral dilution of 
metallic features of old stars caused by young stars, a power-law or a combination of both. 
Indeed, in several of the plots below Wk is used with this function, since, for most practical 
purposes, Wk and xpc are equivalent (Fig. 3). 

4. Emission Lines 

Circum- nuclear starbursts in Seyfert 2's act as a second source of ionizing photons 
besides the hidden AGN, and as such are bound to contribute at some level to their emission 
line spectrum. In this section we examine our sample in search of signs of this contribution, 
analyzing the equivalent widths of strong emission lines (§4.1), emission line ratios (§4.3) and 
line profiles (§4.4). An investigation of line-FC relations is also presented (§4.2 and §4.5). 
When referring to emission line equivalent widths measured with respect to the FC we will 
use the notation W[ c = W^ bs /xfc to distinguish them from the observed ones, denoted by 
W? bs . 

A 

4.1. The Equivalent Widths of H/3, [OIII] and [Oil] 

The effect of a circum-nuclear starburst on the W^ of Seyfert 2's is easily predicted by 
considering the inverse situation: a pure starburst to which we add a Seyfert 2. The latter 
will surely add more to the H(3 flux than to its underlying continuum, since the non-stellar 
FC from the active nucleus is only seen periscopicaly and hence strongly suppressed in the 
scattering process, whereas photons from the Narrow Line Region (NLR) do not suffer such 
suppression because they are not obscured by the pc scale dusty torus which surrounds the 
nucleus. The net effect of this superposition is that Seyfert 2's should have a larger ratio of 
line per FC photons than in pure starbursts. 

This expectation is totally confirmed in the left panels of Fig. 4, where we compare the 
distribution of W^ for our Seyfert 2 sample with that of the 35 Starburst galaxies in the 
Starburst sample. These histograms reveal a clear offset in W^ between Seyfert 2's and 
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Fig. 4. — Top: Distribution of the equivalent widths of H/?, [OIII] and [Oil] for the sources 
in the Seyfert 2 sample. The open histogram represents all galaxies in the sample, whereas 
the shaded regions indicate the composites only. Bottom: As above for galaxies in the 
Starburst sample. The shaded histogram marks galaxies observed through r ap < 1 kpc 
apertures. All W's are measured with respect to the Featureless Continuum, whose strength 
xpc was estimated from our synthesis analysis. In the H/3 and [OIII] plots the FC strength 
is evaluated at 4861 A, whereas for [Oil] we evaluate xfc at 3660 A. 
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stellar powered emission line sources. The median is a factor of 5 larger in the Seyfert 
2 sample than in pure Starbursts (129 compared to 26 A). 

Given their dual nature, we intuitively expect composite systems to exhibit W^ in 
between those of "pure" Seyfert 2's and Starbursts. This expectation is also confirmed. The 
composite Seyfert 2/starburst systems in our sample, marked by the filled region in Fig. 4a, 
tend to populate the low end of the Wf^ distribution in Seyfert 2's, overlapping with the 
high W^p Starbursts (Fig. 4b). 

The same considerations, of course, apply to other emission lines. In fact, given that 
[OIII]/H/3 is much larger in Seyferts than in stellar powered systems, we expect an even 
clearer separation between these two types of objects in terms of W^jjjy This is confirmed 
in Figs. 4c and d. The difference between median W FC, s, which was a factor of 5 for H/3, 
is now 32-fold, with 44 A for Starburst galaxies and 1431 A in the Seyfert 2 sample. As for 
H/3, composites are skewed towards low W^j^s (Fig. 4c). The same effects are identified 
in [Oil] (Figs. 4e and f). 

The large offset in W FC, s between Seyfert 2's and Starburst galaxies is not an artifact 
of aperture differences between the two samples. One does not expect drastic changes in the 
line per continuum photon ratio between nuclear and off-nuclear star-forming regions, so the 
distribution of W FC, s for Starburst galaxies should remain roughly the same for small and 
large apertures. This is demonstrated by the shaded portions of the histograms in Fig. 4b, 
d and f, which mark Starburst galaxies observed through physical apertures r ap < 1 kpc. 
Aperture effects become important for Seyfert 2's, since there the contribution of off-nuclear 
HII regions affects the line/FC proportion, diluting the higher W FC of the nucleus. This 
effect is detected for the objects in common between the Seyfert 2 and NLAGN samples. 
Mrk 477, for instance, has = W^S/xpc — 35/0.47 = 75 A for the observations of 

SBMQ95, whereas with our nuclear spectra we obtain 92/0.87 = 107 A. The same effect 
occurs comparing the nuclear with the "whole aperture" spectra of LK95, which integrate 
over the entire galaxy. 

Since in §3 we found that composites tend to have a strong FC, we may expect them to 
be well separated from "pure" Seyfert 2's in diagrams involving xpc and these emission line 
W FC 's. This is illustrated in Figs. 5a-c, where we plot the W FC 's of H/3, [OIII] and [Oil] 
against Wk, here representing the FC strength. Note that all composites are located to the 
left of the W K = 10 A line indicated in the plots. (A slightly more stringent constraint of 
W K < 8.5 A can be adopted by restraining our definition of composites to systems which 
exhibit signs of young stars, thus excluding the post-starbursts NGC 5643, Mrk 78 and ESO 
362-G8.) The W FC, s by themselves do not segregate composite from "pure" Seyfert 2's as 
efficiently as Wk, but sources with the weakest emission lines are predominantly composites, 
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Fig. 5. — Equivalent widths of H/3 (top), [OIII] (middle) and [Oil] (bottom) with respect to 
the FC against Wk- Sources from the Seyfert 2 sample are in the left (a-c), whereas plots in 
the right (d-f) show the Starburst and Merger samples. Symbols follow the same convention 
as in Fig. 2. The location of Seyfert 2's in the Wk X W^p diagram (a) can be used to assess 
the degree of difficulty to identify starburst features: Sources to the lower left (weak Wk and 
W^p) have a strong, starburst dominated FC, whereas we do not know with certainty what 
dominates the weak FC sources in the top right of the diagram. 
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so the vertical axis provides some extra diagnostic power. The "pure" Seyfert 2 at Wr = 7.4 
A, trespassing the W K < 10 A zone of composites, is NGC 1068, which has W^fi^ = 68, 
WEfijn = 1334 and = 136 A. This prototype Seyfert 2, which shows no signs of recent 
star-formation close to the nucleus, is located either among or close to the composites in 
all diagrams presented in this paper. We will return to this issue in §7. The panels on the 
right side of Fig. 5 show the results for sources in the Starburst and Merger samples for 
comparison. 



4-1.1. The difficulty in identifying circum-nuclear starbursts 

Once again we must draw attention to contrast effects which may lead us to classify 
Seyfert 2's with weak circum-nuclear starbursts as "pure" systems. Part of the difficulty in 
identifying starburst features in Seyfert 2 nuclei with strong emission-lines comes from the 
fact that as Wj^p increases the high order HI Balmer absorption lines of massive stars become 
increasingly filled up by emission. In the absence of far-UV spectra (available for very few 
sources) and the short-lived WR bump, this effect effectively prevents us from seeing the 
starburst in the near-UV to optical range. 

This contrast effect comes on top of the difficulty in identifying circum-nuclear starbursts 
in systems dominated by the old stellar population (§3.2). The combination of these two 
effects may be at least partly responsible for the horizontal and vertical separations between 
composites and "pure" Seyfert 2's seen in Fig. 5a-c. The location of a Seyfert 2 in the Wk- 
Wffp plane can therefore be read as an empirical measure of the "degree of difficulty" in the 
recognition of circum-nuclear starbursts in Seyfert 2's: Sources at the bottom left are easily 
identified as composites, whereas as one progresses to the top right it becomes increasingly 
difficult to discern the starburst features. The usefulness of this concept is illustrated by 
the location of Mrk 3, Mrk 34 and NGC 7212 in Figs. 5a-c, all are around the W K = 10 
A "dividing line" but with W^p and W^ ir ^ larger than for composites. It is no accident, 
therefore, that the nature of their moderately strong FC remains unclear (§3.3.2). 



4.2. The relation between line and FC emission 

It is important to emphasize that the offset in emission line W's between Seyfert 2's 
and Starbursts only appears when these are measured with respect to the FC, as hinted in 
Fig. 5. Our estimate of xfc, which takes into consideration continuum colors and metal 
absorption lines through our semi-empirical synthesis analysis, is thus the key behind this 
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Fig. 6. — Relation between the continuum at 4861 A and H/3 fluxes. The panels differ 
in which component of -F4861 is used in the abscissa, (a) Total observed continuum, (b) 
-^4861 x %old, i-e., the continuum flux due to old stars, (c) F 4861 x xpc, the FC flux. Symbols 
as in Fig. 2. All continuum fluxes are given in ergs -1 cm -2 A -1 . Notice how xfc has the 
property of unveiling a strong line-continuum correlation immersed in the scatter plot on 
the left. 
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significant off-set. 

The role of the synthesis process can be better appreciated in Fig. 6. In Fig. 6a we 
plot F H p against the total observed continuum flux at 4861 A. What is seen is a scatter 
diagram, which is not surprising since the continuum in most sources is dominated by old 
stars and thus has little relation to the line emission. This is confirmed by the fact that the 
scatter increases even further if the continuum fluxes are multiplied by x OLD (Fig. 6b), thus 
isolating the flux due to > 1 Gyr stars. The line versus continuum plot changes dramatically 
when we isolate the FC contribution! Using Fpc — -^Xlli x x fc in the abscissa (Fig. 6c) 
has the remarkable effect of uncovering an underlying order immersed in the scatter plot 
on the left. It is therefore clear that our FC determination method isolated a component of 
the optical continuum which is directly linked to the line emission, lending further credibility 
to the inferred xpc values. The same result is obtained for the sources in our comparison 
samples, i.e., xfc organizes the data along a Fufi-Fpc correlation which is not evident in 
the Fnp-F^-y plane because of the contribution of the old stellar population. 

Possible interpretations of the H/3-FC relation identified in Fig. 6c will be discussed 
in §4.5, after we gather other related results. To substantiate the discussion it is useful 
to compare it to the equivalent relation in Starburst systems and to investigate line-FC 
relations for other transitions. This is done in Figs. 7 and 8, which contain luminosity- 
luminosity versions of Fig. 6c for the Seyfert 2, Starburst and NLAGN samples and for 
R(3, [Oil], [OIII] and Hell. Lines of constant W FC = 10, 100 and 1000 A run diagonally 
across these plots to facilitate their intercomparison. Also to aid the discussion, Table 4 
summarizes the results of a Spearman rank correlation analysis for flux-flux line-continuum 
relations. The values listed are the probabilities of chance correlation (ps), so small values 
indicate significant correlations. 

Significant line-FC correlations are also obtained for [Oil] , [OIII] and Hell in the Seyfert 
2 sample (Figs. 7c, 8a and 8c, Table 4). Note, however, that the scatter in the [OIII] and 
Hell plots is visibly larger than for H/3-FC (Fig. 7a). It is also clear that a substantial 
part of the scatter in these plots is induced by the offset location of the composites towards 
smaller W FC, s. Statistical confirmation of this effect is provided in Table 4, where these 
correlations are examined for the whole sample and for subsets containing only "pure" or 
composite systems. 

Starburst galaxies also follow a tight H/3-FC relation, with Wj^p = 10-100 A (filled 
triangles in Figs. 7b). This is a long known (Terlevich etal. 1991) and well understood 
behavior of stellar powered systems, for which Fpc actually represents the contribution of 
young stars and hence traces the ionizing flux. Current models (e.g., Leitherer etal. 1999) 
predict Wnp between 400 and 10 A for constant star-formation, different IMF's and ages up 
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Fig. 7. — Relation between the H/3 (top) and [Oil] (bottom) luminosities and the FC lu- 
minosity Lpc = L b s x xfc for the Seyfert 2 (left), Starburst and NLAGN samples (right). 
Note that Lpc is evaluated at 4861 A for the H/3 plot and at 3660 A for [Oil]. Diagonal lines 
indicate, from bottom to top, W FC = 10, 100 and 1000 A. Symbols as in Fig. 2. 
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Fig. 8. — Same as Fig. 7, but for the [OIII] and Hell emission lines. No Hell measurements 
are available for the Starburst and NLAGN samples. 
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to 10 8 yr. These are broadly consistent with the values we obtain, but a detailed comparison 
with such evolutionary models will have to await a cross calibration of our semi-empirical 
synthesis decomposition in terms of the theoretical spectra. As for the Seyfert 2 sample, 
the second best line-FC relation for Starburst galaxies is that with the [Oil] line (Fig. 7d). 
This too is a well known relation, which has in fact been used as an alternative estimator 
of star-formation rates in galaxies (Gallagher etal. 1989; Kennicutt 1992; Tresse etal. 1999; 
Jansen, Franx & Fabricant 2001). 

The location of the NLAGN's in these plots (empty triangles) reflects the mixed nature 
of this sample. Some of them behave just like "pure" Seyfert 2's, most notably NGC 3393 and 
NGC 5506, which stand out in Figs. 7b, d and 8b as the sources with highest W FC 's. Most of 
the other NLAGN's are mixed among Starburst galaxies in these plots. This happens due to 
a combination of the aperture effects discussed in §4.1 and the fact that this sample contains 
several composites and some objects which are more consistent with a LINER classification 
(e.g., NGC 1097, NGC 1433). It is thus clear that these galaxies do not constitute a well 
defined comparison sample. For these reasons, we concentrate on the comparison between 
the Seyfert 2 and Starburst sample. 



4-2.1. Extinction 

One would like to believe that extinction does not play a relevant role in these line-FC 
diagrams since both axes correspond to similar wavelengths. In Starburst galaxies, however, 
it is known that the continuum suffers ~ 2 times less extinction than the emission line region 
(Calzetti, Kinney & Storchi-Bergmann 1994). This effect can reshuffle the data points in 
the plots above, potentially changing the overall appearance and strength of the line-FC 
correlations. 

An estimate of the effects of extinction is presented in Fig. 9 for the FL3-FC relation 
using the Ha/H/3 values listed in Table 2. Fig. 9a presents the data with Lpc and Lh@ 
dereddened by the same amount, whereas in 9b we halved the optical depth to the FC. Two 
effects are clear when comparing these plots with their uncorrected counterpart (Fig. 7a): 
First, with the extinction correction, composites are more segregated from "pure" Seyfert 
2's in terms of L H p, with median values which differ by an order of magnitude (median 
L H p = 1.6 x 10 41 and 1.7 x 10 40 ergs -1 for composites and "pure" Seyfert 2's respectively). 
This preference of composites for luminous systems will appear again when we discuss their 
far IR luminosities in §5. Second, a differential line/FC extinction correction only tightens 
the FL3-FC correlation! The effect is such that it practically erases the offset in Wj^p between 
composite and "pure" Seyfert 2's. 
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Fig. 9. — (a) Dereddened FC versus H/3 relation for the Seyfert 2 sample after correcting 
both axis for the extinction implied by Ha/H/?. (b) As above, but assuming the continuum 
is only half as extincted as the emission lines (Calzetti etal. 1994). Symbols and lines as in 
Fig. 7. 
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The caveats in these results are that the Ha/H/3 values used were not obtained from 
the same spectra we analyze (§2) and that we do not know whether the r FC ~ T H p/2 result 
applies to AGN. At any rate, the larger Ha/H/5 of composites, which is responsible for this 
regrouping, is consistent with UV (Heckman et al. 1995) and X-ray (Levenson, Heckman 
& Weaver 2001a) indications that circum- nuclear starbursts in Seyfert 2's are substantially 
reddened, and may significantly contribute to the extinction of the AGN component in 
such composite systems. This "excess" of dust is probably associated with the same gas 
which feeds the starburst activity, but more work will be required to address this issue 
quantitatively. 



4.3. Do starbursts participate in the ionization of the gas in Seyfert 2's? 

The results reported so far show how the presence of circum-nuclear starbursts manifest 
itself in the FC strength and the equivalent widths of emission lines. These results per se 
do not necessarily imply that the starburst contributes to the ionization of the line emitting 
gas in Seyfert 2's, since it may act mostly as a source of FC, diluting an AGN powered H/3 
emission and thus shifting points to the right, but not much upwards in Figs. 6 — 9. One way 
to establish whether the starburst participates in the ionization of the gas is to investigate 
the ratio between HeIIA4686 and H/5. Unlike AGN, starbursts do not produce significant 
radiation above 54 eV, so they show little Hell emission. If starbursts in AGN contribute to 
H/3 this ratio will tend to be smaller in composites than in "pure" Seyfert 2's. 

This is confirmed in Figs. lOa-b, where we plot HeII/H/3 against Wk and xold, both 
empirical indicators of compositeness (§3). Whereas the whole sample spans the 0.1-0.8 
interval in HeII/H/5, all composites have HeII/H/3 < 0.4. Many "pure" Seyfert 2's also have 
Hell this weak, but their stronger Ca II K and hence larger xold clearly separates them 
from composites. Starburst galaxies would be located in the same horizontal range spanned 
by composites in Fig. lOa-b, but down to weaker Hell (a line which is often not detected 
in starbursts). Composites also tend to be less excited in terms of [OIII]/H/5, as shown in 
Figs. lOc-d. The effect of the stellar absorption feature underneath the H/3 emission-line was 
not corrected for in this analysis. Since composites are more likely to have stronger Balmer 
absorption lines than objects dominated by old stars, this correction would only strengthen 
our conclusion, shifting the composites to even smaller HeII/H/9 and [OIII]/H/3 ratios. 

Either line ratio can therefore be used as an auxiliary diagnostic of compositeness. 
This is illustrated by the dotted lines in Fig. 10, which isolate composites from "pure" 
Seyfert 2's almost completely! (Again, the odd ball is NGC 1068.) This shows a welcome 
consistency between diagnostics of compositeness based on continuum and stellar features 
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Fig. 10. — "Excitation" ratios plotted against indicators of circum-nuclear starbursts for 
the Seyfert 2 sample. Filled circles mark composite systems and empty circles mark "pure" 
Seyfert 2's. Dotted lines delimit regions occupied by composites, and may be used as em- 
pirical diagnostics of compositeness. 
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and the expected impact of circum-nuclear starbursts upon the emission line ratios in Seyfert 
2's. Note, however, that while all known composites have low excitation, the converse is 
not true, since many "pure" Seyfert 2's lie in the HeII/H/5 < 0.4 and [OIIIj/H/3 < 15 range 
defined by composites. In terms of emission line ratios, these galaxies could well have circum- 
nuclear starbursts, but which do not make it to the composite class due to the contrast effect 
discussed above. 

It is hard to see how to account for these results in a "pure AGN" scenario. The shape 
of the ionizing spectrum, the proportion of matter to ionization bounded clouds and the 
ionization parameter are just some of the factors which define the excitation level of the NLR, 
and source-to-source variations of these properties can account for a wide spread in HeII/H/5 
and [OIII]/H/3 (e.g. Viegas-Aldrovandi 1988). Yet, none of these intrinsic properties explains 
the "zone of avoidance" for strong FC and high excitation in Fig. 10. If anything, a strong 
FC would suggest a higher ionization parameter and perhaps a harder ionizing spectrum, 
which would produce higher excitation, contrary to what is observed. Overall, the simplest 
interpretation of these data is that part of the emission lines in composites originates from 
photoionization by a circum-nuclear starburst. This contribution is strongest in H[3, since 
starbursts are inefficient producers of Hell and [OIII] compared to AGN. 

It is clear that in a starburst + AGN mixture of emission lines the starburst compo- 
nent can dominate H/3 without moving the source outside of the AGN region in diagnostic 
diagrams (e.g., Hill etal. 2001; Levenson etal. 2001b). For instance, using typical [OIII]/H/5 
ratios of 10-20 for Seyfert 2's and 0.3-1 for metal-rich starbursts, the total ratio still lies 
within the [OIII]/H/3 > 3 domain of Seyferts for as much as 70-90% of H/3 powered by 
the starburst! Therefore our conclusion that photoionization by massive stars provides a 
substantial part of L Hf3 in composites is not in conflict with their Seyfert 2 classification. 
In fact, this conclusion was already implicit in the off-set in between composites and 

Starburst galaxies (Fig. 4). The values in composites cover the range from 30 to 200 

A, which may be represented by its geometric mean of ~80 A. The corresponding value in 
Starburst galaxies is ~30 A, representing the = 10-100 A interval. Since we know that 
the FC in composites is dominated by the starburst component, these values imply a typical 
starburst contribution of 30/80 ~ 40% to H/5. 

These results imply that H/3 is partially powered by the starburst in composite systems. 
Unfortunately, as for Lpc, it is not possible to unambiguously separate the starburst and 
AGN shares of Lh@ without extra information, such as assuming typical line ratios or equiv- 
alent widths. For instance, taking the mean [OIII]/H/3 of 13 for our "pure" Seyfert 2's and 
8 for composites, one obtains a starburst contribution to H/3 of 40% (precisely what we have 
derived above comparing composites and Starburst galaxies in terms of W£g). Illustrative 
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as these global estimates are, it is clear that both composite and "pure" systems present a 
range of properties. For the composites below NGC 5135 ([OIII]/H/5 = 4) in Fig. 10c and d, 
for instance, the same exercise would yield more than 70% of H/3 powered by the starburst. 
Much therefore remains to be learned from a more detailed object by object analysis. 

4.4. Emission Line Profiles: Kinematical separation of the starburst and AGN 

One prospect for disentangling the starburst and AGN contributions to the emission 
lines is to study their profiles in search of differential kinematical signatures of these two 
components. We have thus performed a line profile analysis for the Seyfert 2 sample, but 
only for the Kitt Peak data, which have higher spectral resolution. 

Most of the composites have [OIII] and H/3 profiles which can each be described in terms 
of two components. Interestingly, the narrower component often has a lower excitation 
than the broader one (Table 5), qualitatively consistent with them being powered by a 
starburst and AGN respectively. In fact, the excitation of the broad component (mean 
[OIII] /H/5 = 10.7, standard deviation = 2.8) is very similar to the excitation of "pure" 
Seyfert 2's. However, the excitation of the narrow components ranges between 2.7 and 11.6. 
Meanwhile, the excitation in Mrk 1066, Mrk 1073, NGC 5135, NGC 7130 and IC 3639 is 
similar to the excitation in prototypical nuclear starbursts (e.g. NGC 7714, Gonzalez-Delgado 
etal. 1995), while in Mrk 1, Mrk 463E, Mrk 477 and Mrk 533 the excitation is higher and 
similar to the excitation in young and very highly excited HII regions (e.g. NGC 2363, 
Gonzalez-Delgado etal. 1994). The difference in excitation between these two sub-groups 
within the composites is suggestive of evolutionary effects, where the youngest systems are 
the ones with highest excitation. This is consistent with other properties of composites, 
such as Wfjp , the fraction of total optical light provided by intermediate and young stars 
(xtnt/xfc), and the detection of Wolf-Rayet features (see §7.3). 

If this interpretation is correct, i.e., the narrow and broad emission line components 
are produced by the starburst and AGN respectively, the ratio of the narrow to total line 
flux gives an estimate of the starburst contribution to the ionization of the gas. For H/5, 
this ratio ranges between 37% (in Mrk 463E) and 80% (in Mrk 1073 and NGC 5135), in 
agreement with our coarser estimates in §4.3. These values suggest a significant impact of 
the circum-nuclear starbursts upon the emission line ratios in Seyfert 2s. 
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Table 5. Line Profile Analysis of Composites 



Name ( ™) ( ™) [%] ^ 

V ^ / narrow V ^ / broad F »0 



narrow 



Mrkl 


11.6 


11.5 


38 


38 


Mrk273 


3.1* 








Mrk463E 


8.9 


8.3 


37 


39 


Mrk477 


8.8 


8.7 


58 


59 


Mrk533 


9.4 


12.5 


54 


46 


Mrkl066 


2.9 


5.7 


70 


55 


Mrkl073 


4.6 


15.4 


80 


54 


NGC5135 


2.7 


10.5 


80 


54 


NGC7130 


2.7 


11.7 


60 


25 


IC3639 


5.4 


12.0 


66 


47 



*The nuclear spectrum of Mrk 273 corresponds to the central 2.1x1.5 arc- 
sec around the continuum maximum at 4450 A. This maximum is shifted 
with respect to the maximum of the nebular emission lines. The [OIII] map 
obtained by Colina, Arribas & Borne (1999) shows an off-nucleus Seyfert 2 
nebula 4 arcsec south with a [OIII]/H/3 ~ 8, similar to the excitation of the 
broad component in the other composites. 

Note. - - The principal properties of the two kinematic components fitted 
to the H/3 and [OIII] A5007 emission lines in the composites Seyfert 2 nuclei. 
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4.5. What drives the line-FC correlations in Seyfert 2's? 

The FC strength, W^c's, emission line ratios and profiles all consistently reflect the pres- 
ence of circum- nuclear starbursts in Seyfert 2's and can therefore be used as the empirical 
diagnostics of compositeness which we set out to identify in this paper. The line-FC correla- 
tions presented in Figs. 7 and 8 do not play any direct part in such diagnostics. Nevertheless, 
because of their high statistical significance, potential relevance and controversial history, in 
this section we speculate on what might be driving them. We do this by contrasting two 
extreme views on the origin of the FC in "pure" Seyfert 2's. 



4-5.1. Starburst- dominated FC 

Since we have established that circum-nuclear starbursts make substantial contributions 
both to the FC and H/3 emission, and given that these two quantities are causally linked in 
starburst systems, the presence of circum-nuclear starbursts in Seyfert 2's naturally leads to 
a Hf3-FC correlation. Even if the AGN components of F FC and F H/3 are uncorrected, the 
starburst portions should suffice to maintain a H/5-FC correlation in the combined AGN + 
starburst data. 

One extreme interpretation of Fig. 7a is thus that the H/5-FC correlation is essentially 
driven by the presence of circum-nuclear starbursts. If this is to apply to all sources, the 
FC must be dominated by the starburst even in "pure" Seyfert 2's. Varying proportions 
of starburst to AGN ionizing power can be invoked to account for the few "pure" Seyfert 
2's whose H/3/FC ratios exceed the maximum W^g = 200 A reached by composites. This 
can be explained either by scaling up the nuclear ionizing source, thus increasing Fnp while 
keeping F FC constant, or by scaling down the starburst, which also increases Wj^p since 
F FC would be more affected than F HB . Such adjustments can even be dispensed with if one 
allows for differential extinction to the line and FC, since we have seen in Fig. 9b that this 
alone substantially reduces the W§p differences between composite and "pure" systems. 

But how do the [OIII] (Fig. 8a) and Hell-FC (Fig. 8c) correlations fit into this model? 
Table 4 shows that these correlations, while poorer than the one between F H3 and F FC , are 
still significant at the > 99% confidence level, and that the reason for the larger scatter is 
the small WEfjjn and W^n °f some 5 or 6 composites. Such a scatter is in fact expected in 
a starburst dominated FC model, since the [OIII] and (especially) the Hell line fluxes trace 
the AGN power, whereas Fpc is a measure of the starburst power, so that and W^n 

are indicators of the varying contrast between these components. Decreasing the power of 
the starbursts in composite systems by factors of 2-3 while keeping the AGN constant would 
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move composites along ~ horizontal lines to the left in Figs. 8a and c, placing them among 
"pure" Seyfert 2's. 

While it is easy to see why composites deviate in W^ HI ^ and W^ H) it is not clear why 
"pure" Seyfert 2's should define line-FC correlations as good as those revealed in Figs. 7 
and 8 (see also Table 4). If their FC is indeed dominated by a weak starburst, then the 
fact that their AGN-dominated Fyoni] and Fneii scale with Fpc would indirectly imply the 
existence of a proportionality between the starburst and AGN powers, i.e., that more powerful 
starbursts occur in more powerful AGN! Further evidence for this scaling is presented in §5, 
§7.2 and GD01. 



4.5.2. AGN- dominated FC 

An opposite model would be to postulate that the FC is dominated by starbursts only in 
composites, while in "pure" Seyfert 2's the FC is dominated by scattered light, so these two 
classes should not be mixed when discussing line-FC relations. This interpretation would 
lead us to the intriguing conclusion that the scattering efficiency, e, which links the observed 
FC1 to the intrinsic nuclear FC via Lpci — tLpco, does not vary substantially among Seyfert 
2's. The surprise comes from the fact that e depends on the geometry of the mirror and the 
optical depth of scattering particles (Miller, Goodrich & Mathews 1991; Cid Fernandes & 
Terlevich 1995). Since there is no a priori reason why e should be similar for all Seyfert 2's, 
we would expect a wide spread of W FC 's, contrary to the relatively narrow ranges observed, 
particularly if composites are excluded. 

This argument has in fact been pursued by Mulchaey et al. (1994) to explain the absence 
of line-FC correlations in their data, from which they claim a consistency with the unified 
model. As acknowledged by Mulchaey etal., however, it is possible that, just as we find 
here, a substantial part of the scatter in their line-FC plots is associated with contamination 
by starburst activity. That this is more than a possibility is illustrated by the fact that the 
Seyfert 2 list of Mulchaey etal. (1994) includes 12 of our 15 composites! It also includes 
other prime composite suspects, most notably NGC 6221, a source for which Levenson etal. 
(2001b) find a complete dominance of the starburst, with the AGN shining through only 
in hard X-rays. As a whole, their sample is similar to the NLAGN sample of SBMQ95 
in terms of the mixed types of sources. Indeed, as seen in Figs. 7b, d and 8b, no line- 
FC correlations are present for this heterogeneous sample. Furthermore, contamination by 
starbursts is certainly larger in the Mulchaey et al. data than in ours, since they trace the FC 
by the far-UV flux collected through the large aperture of the IUE spectrometer, a procedure 
which follows the premise that starlight contamination is irrelevant in the UV. This premise 
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was proven wrong by Heckman etal. (1995), who have demonstrated that, on the contrary, 
scattered light from a hidden Seyfert 1 is a minor contributor to the IUE spectra of Seyfert 
2's, which are more typical of reddened starbursts. 

At least in relative terms, it is clear that our observations are better suited to investigate 
line-FC relations. Our results essentially recover the classic Yee (1980) and Shuder (1981) 
correlation between line and FC luminosity, which extends over more than 5 decades in 
luminosity from Seyfert 2's to QSOs. Looking at their original plots one sees that Seyfert 
2's show a larger scatter around the correlations than type 1 sources, but they still broadly 
follow the correlations. While it is true that the interpretation of their line-FC correlations 
was left orphan with the advent of the unified model, their reality has not been convincingly 
challenged as yet. Similarly, barring unidentified selection effects, it is hard to deny the 
reality of the line-FC correlations seen in Figs. 7 and 8. 

In a Lpc = Lpci model equals W^^e -1 , where W^ 1 is the equivalent width 

which would be measured from non-obscured lines of sight to the Seyfert 1 nucleus. Ob- 
servationally, the narrow H/5 in Seyfert l's roughly spans the range between W^ 1 ~ 5 to 
75 A (Goodrich 1989; Rodriguez- Ardila, Pastoriza & Donzelli 2000). This spread is already 
comparable to the 70 < < 560 A interval span by "pure" Seyfert 2's in our sample, 

which does not leave much room for a spread in e. Combining these intervals we obtain an 
acceptable range for e of roughly 0.01-0.1. (As a reference, for NGC 1068, the only source 
for which the scattering has been modeled in detail, e ~ 0.015; Miller etal. 1991). A "well 
behaved" e is thus required in order to explain the H/5-FC correlation for "pure" Seyfert 2's. 
This interpretation is analogous to, and just as puzzling as, the interpretation of the con- 
stancy of the W of broad H/3 among type 1 Seyferts and QSOs as due to a roughly constant 
covering fraction of broad line clouds in all sources (Yee 1980; Shuder 1981; Osterbrock 1989; 
Binette etal. 1993). While viable, this model obviously requires some fine tuning. 

Yet another possibility is that the FC carries a sizable nebular component in "pure" 
Seyfert 2's, which would naturally produce line-FC correlations. A purely nebular FC, 
however, would produce of order 2000 A for electron temperatures between 10000 and 
20000 K, much larger than the values we find. A hotter plasma, however, can enhance the 
continuum emission and lower the W FC 's. This is the essence of Tran's (1995c) model for 
FC2: emission from the ~ 10 5 K electron-scattering mirror. If this emission is somehow tied 
to the line emission — which is reasonable since both arise in circum-nuclear gas heated by 
the active nucleus — then a line-FC correlation follows. 

While our personal prejudices tend to favor a starburst-dominated FC interpretation, 
we are not in a position to rule out a role for a non-stellar FC for "pure" Seyfert 2's. Indeed, 
the fact that NGC 1068 mingles so well among composite starburst /Seyfert 2's is a clear 
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demonstration that our analysis cannot resolve the ambiguous nature of the FC in "pure" 
Seyfert 2's. Since we know that several composites have a scattered component (Table 2), the 
correct interpretation of the line-FC relation will certainly involve a mixture of the extreme 
possibilities discussed above. 

5. Far IR luminosities 

In Figs. 7-9 we have seen that, besides having smaller emission line W FC, s, the ten- 
dency of composites to live preferentially among luminous sources allows them to be more 
clearly separated from "pure" Seyfert2's in line-FC diagrams. To evaluate this tendency in 
a reddening independent fashion, Fig. 11a shows the far IR luminosity, computed with the 
60 and 100 /im IRAS bands (cf. Sanders & Mirabel 1996), against Wk for the Seyfert 2 sam- 
ple. The tendency identified by GD01 for composites to be powerful IR emitters is clearly 
seen in this figure. Except for Mrk 1210, ESO 362-G8 and NGC 5643, all composites have 
L FIR > 10 10 L . The difference in median L FIR between composites and "pure" Seyfert 2's 
is a factor of 5. 

Comparison with Fig. lib shows that composites are also more powerful far IR sources 
than typical Starburst galaxies, being more comparable to the merger systems studied by 
LK95. Normal galaxies with little star-formation would populated the large W K and low 
Lfir portion of the plot (Telesco 1988; Soifer, Houck & Neugbauer 1987). 

IR luminous galaxies have a well know preference to be interacting systems (Sanders & 
Mirabel 1996). Fig. 11a therefore suggests that interactions and the associated gas fueling of 
the nuclear regions may be a key ingredient in triggering circum-nuclear starbursts in AGN. 
Indeed, a substantial fraction of our composites (9 of 15) are associated with interacting 
systems or groups (GD01; Levenson etal. 2001a; SB00): NGC 7130, Mrk 1, Mrk 273, Mrk 
463E, Mrk 477, Mrk 533 and IC 3639 are all in interacting systems, while NGC 5135 and 
NGC 7582 belong to groups. The analogous fraction for the the "pure" Seyfert 2's is only 
4 out of 20 (Mrk 348, Mrk 607, NGC 5929 and NGC 7212). This may explain why the 
composites resemble the mergers in Fig. 11. 

Since Wk works essentially as a measure of the starburst strength, Fig. 11 indicates a 
link between the far IR and the amount of recent star-formation in Seyfert 2's. This link is 
investigated in Fig. 12. The abscissa in Fig. 12a is the total observed luminosity at 4861 A, 
whereas in Fig. 12b we filter out the light from old stars by using Lsb = -^4861 x ( x fc + %int)- 
The post-starburst population, which was not included in the line-FC analysis because it does 
not ionize the gas, is now included because it can heat the dust. Furthermore, Table 4 shows 
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Fig. 11. — Far IR luminosity against Wk for (a) the Seyfert 2 sample, and (b) the Starburst 
and Merger samples. Symbols as in Fig. 5. Starbursts with Lpm < few xlO 9 L are generally 
blue compact dwarfs or related objects. Normal galaxies would populate the bottom part of 
the figure (i.e., low L FIR ), with Wk depending on the amount of star-formation. Seyfert 2's 
lie along the upper envelope, with composites skewed to high Lfir- 
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Fig. 12. — Top: Relation between the far-IR and the continuum luminosity at 4861 A for 
the Seyfert 2 sample. In (a) the total L4861 is used in the abscissa, whereas only the FC 
+ post-starburst component is used in (b), and in (c) L 4861 x (xpc + x int) is corrected 
for extinction (using Ha/H/5). Composites are marked by filled circles. Bottom: As above, 
but for the Starburst (filled triangles) and Merger (squares) samples. Abscissa in units of 
ergs^ 1 A -1 . The dotted lines correspond to Lpm/ (vL^^i ratios of 1000, 100, 10 and 1. 
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that the far IR flux is somewhat better correlated with F$b than with Fpc alone. Comparison 
of Figs. 12a and b show how xfc + x int has the property of unveiling a correlation immersed 
in the scatter when the whole continuum is used. An estimate of the effects of extinction is 
presented in Fig. 12c, where we deredden L$b using the Ha/H/3 values listed in Table 2. The 
bottom panels in Fig. 12 show the corresponding results for Starburst galaxies and mergers. 
Comparison of the top and bottom panels shows that Seyfert 2's have Ffir/Fsb ratios larger 
than galaxies in the Starburst sample, but compatible with the LK95 mergers. Despite 
the mixed emission line properties of the sources in the Merger sample, it is reasonable 
to attribute most of their FIR and FC emission to star-formation in a dusty environment 
(Lutz etal. 1999; 1998). This similarity is thus suggestive of the presence of powerful dusty 
circum-nuclear starbursts in Seyfert 2's, consistent with our knowledge of composites. A 
usual caveat when dealing with IRAS fluxes is that one has to allow for the fact that a 
substantial fraction of the far IR light originates well outside the nuclear region. Aperture 
effects may therefore be partially responsible for the larger Ffir/Fsb found for objects in 
the Seyfert 2 than in the Starburst sample. This effect was detected for sources in common 
between the Seyfert 2 and NLAGN samples, the latter having systematically larger F$b 
and thus smaller F fir /F S b- At any rate, it is clear that circum-nuclear starbursts can be 
responsible for the bulk of the far IR emission in Seyfert 2's. This was in fact confirmed in 
the few cases where this issue has been examined more closely (e.g., Mrk 477 Heckman et al. 
1997; see also Rigopoulou etal. 1999). 



6. Near UV surface brightness 

We have shown that the near-UV nuclear light in the Seyfert 2 composites is dominated 
by young and intermediate-age stars. In order to make some causal connection between these 
young stars and the Seyfert phenomenon, we must establish that the properties of this stellar 
population are unusual compared to the nuclei of normal galaxies. Young and intermediate- 
age stars dominate the near-UV continuum in the nuclei of normal galaxies of Hubble type 
Sbc and later (e.g. Heckman 1980). Seyfert galaxies are generally classified as early type 
disk galaxies (e.g. Nelson & Whittle 1996; Ho, Filippenko & Sargent 1997), whose nuclei are 
dominated by old stars. However, the Hubble types of Seyfert galaxies based on ground-based 
images are quite uncertain, since they are generally much more distant than the galaxies used 
to define the Hubble sequence. Indeed, the HST imaging survey by Malkan, Gorjian & Tarn 
(1998) reveals a significant number of Seyfert 2's (including some of our composites) in 
which the HST-based Hubble type is much later than that assigned from earlier ground- 
based images (Storchi-Bergmann etal. 2001). Could it be that the Seyfert 2 composites 
are simply otherwise normal late-type galaxies containing an AGN? The abnormally large 
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Fig. 13. — (a) Distribution of U-band surface brightness for the Seyfert 2 sample, with 
composites indicated by the filled areas, (b) As above, but for normal galaxies from the 
Jansen et al. (2000) catalog, extracted through 250 or 500 pc aperture radii, (c) As above 
but for the Starburst and Merger samples, excluding sources observed through apertures 
larger than r ap = 1 kpc. For the Seyfert 2, Starburst and Merger samples fj,u was measured 
from the F 3660 flux divided by the aperture area. For the Jansen et al. data jiu comes from 
aperture photometry in U. 
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Fig. 14. — Bottom: U-band surface brightness against Wk for the Seyfert 2 sample. Top: 
Same as above after correcting fijj by the extinction using Ha/H/3. Considering that the 
extinction in the continuum can be smaller than that of the emission lines (Calzetti et al. 
1994), panels a and b bracket the correct values for 
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far-IR luminosities of the composites argues against this. In this section we also show that 
the observed UV surface-brightnesses of the Seyfert 2 composite nuclei are much larger than 
the nuclei of normal galaxies (of any Hubble type), but rather similar to typical starbursts. 

We have compared the surface brightnesses in the near UV (~ 3660 A) of all our 
Seyfert 2's with those of normal galactic nuclei based on the imaging data published by 
Jansen etal. (2000). For the normal nuclei (spanning the full range of Hubble types) we 
have computed the mean surface brightness through two metric apertures with radii of 250 
and 500 pc, chosen to represent typical values of the spectroscopic apertures used for the 
Seyfert 2 nuclei (see Fig. 1). This was done converting the U-band surface photometry of 
Jansen et al. (2000) to monochromatic surface brightnesses at 3650 A, including a correction 
for Galactic extinction (from NED). We have also measured near-UV surface brightnesses 
from the spectra of the Starburst and Merger samples. In these samples we have restricted 
our comparison to galaxies in which the projected aperture size of the spectra had a radius 
less than 1 kpc (see Fig. 1). Thus, in all cases we are comparing UV surface-brightnesses on 
similar physical scales (radii of several hundred pc). 

The results are shown in Fig. 13, from which we draw the following conclusions. First, 
the Seyfert 2 composites have much higher central UV surface brightnesses than normal 
galaxies of any Hubble type (by roughly an order-of- magnitude on average). A correction 
for dust extinction would undoubtedly increase the disparity still further. Second, the Seyfert 
2 composites have rather similar central UV surface-brightnesses to the members of the Star- 
burst and Merger samples. When the UV surface-brightness is corrected for extinction (see 
Fig. 14), the Seyfert 2 composites reach values similar to those which would be exhibited by 
the "maximal" starbursts studied by Meurer etal. (1997). Third, the Seyfert 2 composites 
have tend to have higher central UV surface-brightness than the "pure" Seyfert 2's (by an 
average factor of about 2). As we show in Fig. 14, the difference between the two classes is 
larger (factor of ~ 4) when an extinction correction is made. Fourth, the "pure" Seyfert 2's 
overlap with the normal galactic nuclei having the highest UV surface brightnesses. These 
latter tend to be the nuclei of normal early-type (Sb and earlier) galaxies whose high UV 
surface brightness is due to their prominent bulge component. This is consistent with "pure" 
Seyfert 2's being otherwise normal early- type galaxies with only a modest contribution of 
light from an AGN or young stars (as shown by our spectral synthesis). Note that the out- 
lying Seyfert 2 with extremely high surface brightness is the nucleus of NGC 1068 (due to 
the small physical size of the spectroscopic aperture and the strong contribution by scat- 
tered AGN light). Finally, Fig. 14 also shows the expected trend between increasing UV 
surface brightness and an increasing relative contribution to the near-UV light by young and 
intermediate-age stars (as parameterized by Wk)- 
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7. Discussion 

7.1. Diagnostic of compositeness 

Our experiments have consistently shown that Seyfert 2's with circum-nuclear starbursts 
detected by means of either far-UV stellar wind lines, high order Balmer absorption lines or 
the WR bump also present distinctive optical/near-UV continuum, optical emission line, far 
IR continuum, and near-UV surface brightness properties, all in agreement with the expected 
impact of starbursts upon such observables. Given this remarkable consistency, one can use 
the results of our analysis as alternative tools to identify composite sources. 

The following empirical criteria can be used to characterize a Seyfert 2 nucleus as com- 
posite: 

1. Weak Ca II K: W K < 10 A. Or, equivalently, a small contribution of old stars: x ld < 
75%. 

2. Weak emission lines with respect to the FC: < 200 A, Wf c HI] < 2500 A, Wf Q c H] < 
600 A. 

3. Low excitation: HeII/H/3 < 0.4, [OIII]/H/3 < 15. 

4. Large far-IR luminosity: L FIR > 10 10 L . 

5. Composite line profiles, with a narrow component typical of HII- nuclei and a broader 
one with NLR-like excitation. 

6. High near-UV surface-brightnesses on circum-nuclear scales: ~ 10 3 L pc~ 2 uncor- 
rected for extinction and ~ 10 4 L Q pc~ 2 after extinction correction. 

These criteria are laid roughly in order of importance, the Ca II K (or xold) limit being 
by far the most useful. The upper limits for the W FC, s and emission line-ratios quoted above 
were chosen as to encompass all composites in our sample, most of which have substantially 
smaller values than those listed above. As argued above, one cannot base a classification 
on these properties alone, but they are useful to confirm the compositeness for sources 
satisfying the first criterion. Taken together, these empirical rules separate composite from 
"pure" systems very efficiently. 

Naturally, there are border line cases. The moderate FC (15% < x FC < 30%) galaxies 
Mrk 3, Mrk 34 and NGC 7212, which we discussed in §3.3.2, lie either within or close to the 
limits outlined above. The clearest distinction between these galaxies and our composites is 
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in terms of their emission line W's, = 235-321 and Wf Q c in] = 2667-3565 A, typical of 

"pure" Seyfert 2's but larger than in any of our composites, and well away from Starburst 
galaxies. Note, however, that these high W FC, s, coupled with the dominant contribution 
of old stars (xold ~ 77%) are precisely the conditions under which the identification of 
starburst features are most difficult (§4.1.1), so we cannot rule out the possibility that their 
FC is dominated by relatively weak circum-nuclear starbursts. The other moderate FC source 
in our sample is NGC 1068, for which there is overwhelming evidence that the nuclear optical 
FC is predominantly scattered AGN light. The similarities between this prototype Seyfert 
2 and our composites, whose FC are dominated by starbursts, is a reminder of how similar 
starbursts and AGN can be, and a demonstration of the limitations of our diagnostics for 
border line cases. 

Overall, the empirical criteria listed above fulfill our objective of delineating a region 
in a space of easily obtained observables which can be used to indirectly infer the presence 
of circum-nuclear starbursts in Seyfert 2's. We therefore anticipate that their application to 
large samples of near-UV spectra will reveal many more composite systems. In fact, since 
these criteria mix several kinds of properties, it is not at all obvious whether they will hold 
for larger samples. Future applications of these diagnostics will therefore not only serve as a 
means of identifying composite candidates but to verify the consistency of the scheme as a 
whole. It will also be interesting to incorporate other multi-wavelength data in order to help 
disentagling the starburst and AGN components. Valuable tools in this sense were developed 
for the IR part of the spectrum by Miley, Neugebauer & Soifer (1985), Lutz etal. (1998), 
Genzel & Cesarsky (2000), Laurent et al. (2000) and Hill et al. (1999, 2001), among others. In 
the radio part of the spectrum, some of the diagnostics and techniques presented by Kewley 
etal. (2000), Smith etal. (1998a,b), Lonsdale etal. (1993) and Condon etal. (1991) can be 
used, while X-rays signatures of circumnuclear starbursts have been recently investigated by 
Levenson et al. (2001a). 



7.2. A "Starburst oc AGN" connection? 

We have shown that starburst /Seyfert 2 composites have large FL3, FC and far IR 
luminosities within our sample. These are all quantities which may carry a significant, 
even dominant contribution from the starburst itself. In order to investigate whether these 
powerful starbursts are associated with correspondingly powerful AGN we compared the 
[OIII], Hell and 12 /im luminosities of composites and pure Seyfert 2's. These properties 
are less affected by the starburst component and thus provide an empirical measure of the 
AGN intensity. 
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We find that composites are also skewed towards large luminosities in all these prop- 
erties. For instance, 53% of the sources above but just 31% of those below the median 
reddening-corrected L^ ni] of 4.8 x 10 41 ergs -1 are composites. Nearly identical numbers 
apply to Lff e ii, while for L 12/jm , which arguably measures the hot dust heated by the AGN 
(Pier & Krolik 1992; Spinoglio & Malkan 1989), the composite percentages are 65 and 25% 
above and below the median respectively. These asymmetries are comparable to those in 
starburst-sensitive properties. For L H/3 , for instance, we find 65% of composites above the 
median 6.0 x 10 40 ergs -1 , but just 19% bellow it. The corresponding numbers for L FIR are 
71 and 19%. 

Circum- nuclear starbursts in Seyfert 2's, at least those we can detect, are therefore pow- 
erful (bolometric luminosities of > 10 10 L , judging by L FIR ), and surround correspondingly 
powerful AGN. This luminosity link between AGN and circum-nuclear starburst activity can 
be interpreted in several ways. Perhaps the starburst activity stirs up the ISM in a way which 
induces a large rate of gas accretion by the central super-massive black-hole (Colina & Wada 
2000), while less powerful AGN can be feed with more modest accretion rates, dispensing 
the helping hand of a circum-nuclear starburst. While feasible, this does not seem to be the 
only mechanism to produce powerful AGN, since there are also some high luminosity "pure" 
Seyfert 2's in our sample. Or perhaps the lL L starburst oc Lagn" relation extends all over the 
Lagn spectrum, but is only detected for the most luminous AGN because of our difficulties 
to identify weak circum-nuclear starbursts. 

Whatever the correct interpretation, our results strongly suggest that powerful Seyfert 
2's have at least a 1 in 2 chance of harboring detectable circum-nuclear starbursts. This 
prediction can be readily tested by applying the diagnostics developed in this and our pre- 
vious work to near-UV spectra of such systems, some of which have luminosities rivaling 
those of quasars (de Grijp etal. 1992). Extending the starburst- AGN connection studies 
towards the low luminosity end, on the other hand, will require alternative techniques. High 
spatial resolution (i.e., HST) near-UV spectroscopy offers the best hope of maximizing the 
contribution of circum-nuclear star-forming regions with respect to the old stellar population 
background. This technique was recently applied to the LINER NGC 3507 by Shields et al. 
(2000), where young stars previously undetected from the ground were clearly revealed with 
HST. 



7.3. The meaning of compositeness 

The composite vs. "pure" Seyfert 2 classification which runs throughout this paper 
does not necessarily reflect an intrinsic dichotomy among type 2 Seyferts. In the spirit of 
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Fig. 15. — (a) "Evolution" of the H/5 equivalent width with respect to the continuum flux 
due to young and intermediate age stars for the composites in the Seyfert 2 sample. The 
Xint/xfc ax is is an empirical age indicator, which roughly maps onto a to a few xlO 8 yr 
age interval within the limits above, (b) and (c) As above, but for the AGN sensitive Hell 
and [OIII] lines. The arrow at the bottom right of (a) and (b) marks the horizontal position 
of ESO 362-G8, for which we did not measure H/3 nor Hell emission. 
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unification, one would rather prefer to look at these two classes as representing extremes of 
a continuous distribution of some underlying physical parameter. This view is warranted by 
the continuity and overlap in observed properties between composites and "pure" Seyfert 2's 
found in this and other studies (Levenson etal. 2001a; Veilleux 2001 and references therein). 

At least two ingredients must be involved in any physical description of a starburst- 
AGN connection, one of which is some measure of the relative starburst to AGN intensities, 
such as the ratio of star- format ion to accretion rate, or the black-hole to stellar mass ratio. 
Of the properties investigated in this paper the W s, mainly of [OIII] and Hell, and the 
excitation ratios [OIII]/H/9 and HeII/H/5 are in principle the ones more closely associated 
with this "starburst-to-AGN ratio" . The fact that as a whole "pure" Seyfert 2's behave more 
like AGN in terms of these indicators than composites is thus consistent with them having 
a smaller starburst-to-AGN ratio. The strongest contrast effect, however, is between the 
starburst and old stellar population, as deduced from the fact that all composites but only 
one "pure" Seyfert 2 have xold < 75%, a limit which can be expressed in a more empirical 
way by Wk < 10 A. The Ca II K line originates in bulge stars, being neither a tracer of 
starburst nor AGN activity, so this threshold is clearly an observational limitation which 
prevents us from identifying circum-nuclear starbursts not powerful enough to dilute Wk 
below 10 A. The "pure" category ought to contain at least some such systems, and in fact 
several "pure" Seyfert 2's in our sample, particularly the less luminous ones (see Fig. 7), 
have W FCl s and excitation ratios within the range spanned by composites. 

The other ingredient is evolution. The combination of starburst-to-AGN ratios and 
evolution can in principle account for the whole range between composite and "pure" systems. 
While our coarse, 3-component population synthesis analysis does not provide a detailed 
description of the evolutionary history of circum-nuclear starbursts in Seyfert 2's, in §3.2 
we have seen that the relative strengths of xpc an d %int in composites bear a very good 
correspondence with the evolutionary status deduced from a much more detailed spectral 
analysis (SB00; GD01). The ratio xint/xfc can thus be used as an empirical measure of 
the time elapsed since the onset of star-formation, allowing us to broadly assess evolutionary 
trends in the data. 

Fig. 15 shows that younger systems have stronger H/3, Hell and [OIII] emission lines, 
which we tentatively interpret as evidence for evolutionary effects. Only composites are 
plotted, since only for them we can safely interpret the FC component as due to young 
stars and hence read xfc/ x int a s an age indicator. Also, note that W's in this figure are 
measured with respect to the whole "starburst flux" F SB = F obs (xpc + %int)- This choice of 
normalization facilitates the comparison with spectral evolution models such as Starburst99 
(Leitherer etal. 1999), particularly for systems undergoing star- formation over a period of 
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order 10 s yr. Since H/3 contains a strong contribution from gas ionized by the starburst, 
Fig. 15a roughly traces the evolution of the starburst. The Hell and [OIII] to F SB ratios, 
on the other hand, are measures of the relative power of the AGN and the starburst, and so 
Fig. 15b and c can be used to assess the evolution of the AGN with respect to that of the 
starburst. 

As two extremes, we consider both an instantaneous burst and constant star-formation 
to describe the starburst. In the former case, the optical continuum due to the young stars 
peaks in luminosity about 3 Myr after the burst, and then fades by a factor of 3 by 10 Myr 
and by another factor of 10 by 50 Myr. The systematic decrease in Wueii and W[oiii\ in the 
more evolved starbursts would therefore require that the AGN fades more rapidly than this. 
In this instantaneous burst scenario, the AGN lifetime must be very short in the composite 
Seyfert 2's (< few xlO 7 years). Since the AGN should disappear before the intermediate-age 
stars, the descendants of the Seyfert 2 nuclei would be part of the post-starburst population 
(the "E+A" galaxies). Thus, the only plausible evolutionary link between the composite and 
"pure" Seyfert 2's would seemingly require that latter are the evolutionary precursors to the 
former. This does not seem physically appealing, among other reasons, because the AGN's 
in the "pure" Seyfert 2's are significantly less powerful on-average than in their supposed 
descendants (the youngest composites). Moreover, an instantaneous burst description is 
clearly not appropriate for sources like NGC 5135 and NGC 7130, which contain a mix of 
young and intermediate age stars. 

Alternatively, we can consider a more long-lived (3> 10 7 years) event in which the star- 
formation remains roughly constant. The Wnp values predicted by Starburst 99 for constant 
star formation vary from 150-450 to 10-60 A between and 10 8 yr for different initial mass 
functions and solar metallicity. This range covers the values in Fig. 15a. The ~ 10-fold 
decrease in Wnp during this period is also in agreement with the data. In this model, the 
luminosity of the optical continuum steadily grows with time (e.g. by a factor of 3 to 4 
between 10 and 300 Myr after the start of the burst). The systematic decrease in Wneii 
and Woiii in the more evolved starbursts would therefore imply that the AGN's luminosity 
remains relatively constant over this time-scale. Long time-scale bursts are indeed seen in 
pure starburst galaxies, where ground based spectra usually integrate over regions containing 
several young, but non-coeval, super- associations (Meurer et al. 1995; Gonzalez Delgado et al. 
1999; Langon etal. 2001). Individual starburst knots have been identified in some of our 
composites (Gonzalez Delgado etal. 1998), so we expect this averaging effect to happen in 
Seyfert 2's too. In this scenario, both the starburst and AGN are relatively long-lived, and 
the AGN could well outlive the starburst. As the star- formation ceases, the starburst will 
fade first in H/3, then in the optical continuum, eventually reaching the point where only 
scattered FC and AGN-powered lines survive. With the drop in the starburst-continuum 
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the FC becomes weaker and W K less diluted. Galaxies like ESO 362-G8, NGC 5643 and 
Mrk 78, with their weak FC but pronounced post-starburst populations may represent this 
intermediate stage. [OIII] and Hell will be much less affected than B.f3 by the death of 
massive stars, so excitation ratios will go up. Starburst-heated dust emission will subside, 
decreasing Lpm- All these suggest that composites will end up as "pure" Seyfert 2's. 

Clearly, the above is quite speculative. Eventually, it may be possible to test these 
ideas empirically with improved age dating techniques which will allow us to evolve detected 
starbursts back and forth in time and thus to evaluate whether the age of circum-nuclear 
starbursts may be used as a clock for nuclear activity. As in any other astrophysical context 
where evolution is important, statistically significant samples will be necessary to trace all 
relevant stages. 

8. Conclusions 

We have examined a sample of 35 galaxies previously studied in a series of papers 
which addressed the occurrence of circum-nuclear starbursts in Seyfert 2's and implications 
for the "starburst-AGN connection". The detailed spectral analysis performed in these 
previous papers was used to classify the nuclei into 15 composite and 20 "pure" Seyfert 
2's, with composites defined as nuclei in which large populations of stars with ages of 10 6 - 
10 8 yr have been conclusively detected either through UV stellar winds features, high order 
Balmer absorption lines in the near UV or the WR bump. Using this classification as a 
guide, we investigated continuum colors, absorption line equivalent widths, emission line 
properties, far-IR luminosities and near-UV surface brightnesses for the sample galaxies 
and comparison samples of normal and starburst galaxies, with the primary goal of verifying 
whether these properties reflect the presence or otherwise of circum-nuclear starbursts. Since 
this consistency was verified, these properties, all of which are more easily measured than 
those originally used to detect these starbursts, offer alternative empirical diagnostics of 
compositeness which can be applied to larger and more distant samples. 

Our main results can be summarized as follows. 

(1) We have applied a semi-empirical population synthesis method which effectively pro- 
vides a bi-parametric description of each galaxy in terms of old (> 10 9 yr), post- 
starburst (10 8 yr) populations and a "Featureless Continuum" which represent the 
combined contribution of < 10 7 yr stars and a non-stellar FC. In practice, this represen- 
tation requires little more than the equivalent width of the Ca II K line, and empirical 
calibrations were presented which allow an estimation of the (xold,%int,xfc) vector 
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from Wr and the -F 3660 /-F402o color, both easy to measure observables. This very sim- 
ple scheme proved instrumental in unveiling a number of properties and gives results 
which are consistent with a more detailed spectral analysis of the stellar population 
mixture. 

(2) All starburst/Seyfert 2 composites but only one "pure" Seyfert 2 (NGC 1068) have 
xold < 75% at 4861 A. This threshold can be expressed in a more empirical way 
by Wk < 10 A. Given that our operational definition of composite/ "pure" systems is 
based on the detection or otherwise of stellar features from massive stars, the most 
natural interpretation for this result is that it is due to a contrast effect: Seyfert 2's 
whose circum-nuclear starbursts are overwhelmed by the light from the old stars in the 
host galaxy are not recognized as composite systems. At least some "pure" Seyfert 2's 
must fit this description, i.e., composites with weak starbursts. This may in fact apply 
to most, if not all, "pure" Seyfert 2's, since we cannot determine the nature of their 
weak UV-excess with the data discussed here. 

(3) A related result is that strong FC sources, those with xfc ^ 30%, are invariably 
associated with certified starburst/Seyfert 2 composites. We have argued that this is 
expected, since if scattered light was dominant in such sources then reflected broad 
lines would render them classified as type 1 Seyferts. 

(4) The equivalent widths of emission lines (H/3, [OIII], [Oil]), when measured with respect 
to the FC, are larger in "pure" Seyfert 2's than in Starburst galaxies. Composites have 
W FC, s intermediate between these two extremes. This fits with the idea that circum- 
nuclear starbursts in Seyfert 2's contribute proportionately more to the FC than to 
the line flux, since, while the AGN-powered line emission is seen directly, only a small 
fraction of the corresponding non-stellar FC reaches the observer after scattering. 

(5) Part of the separation in between "pure" Seyfert 2's and composites may be due 
to another contrast effect, in which the high order Balmer lines, the main diagnostic 
of massive stars in the optical, are diluted and eventually filled by nebular emission in 
strong lined sources, preventing them from being recognized as composites. 

(6) Composite systems tend to be less excited in terms of the HeII/H/9 and [OIII]/H/3 
ratios. This indicates that the circum-nuclear starbursts have a direct impact on the 
emission line properties of Seyfert 2's, with typically ~ 50% of H/5 being powered by 
OB stars. 

(7) The dual nature of the ionizing source in composites also leaves a kinematical imprint 
on the emission line profiles, which can be described in terms of a broad component with 
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NLR-like line ratios and a narrower and (usually) less excited one. These components 
are presumably associated with the gas ionized by the AGN and starburst respectively. 

(8) Composite galaxies are luminous in the far-IR, with L FIR ~ 10 10 to 10 12 L and are 
more luminous than "pure" Seyfert 2's and normal galaxies. As a whole, composites are 
rather similar to merger systems in terms of their far-IR properties. About 50% of the 
composites, but only 20% of "pure" Seyfert 2's are in interacting systems, suggesting 
that galaxy interactions spur circum-nuclear starbursts. 

(9) Composite galaxies have near-UV surface-brightnesses in their central-most few hun- 
dred pc that are on-average an order-of-magnitude higher than in normal galaxies of 
any Hubble type. These high surface brightnesses are similar to those in pure starburst 
galaxies. They are also several times higher on-average than in the "pure" Seyfert 2's 
(whose central surface brightnesses in most cases overlap those of normal early-type 
galaxies) . 

(10) These properties were condensed into a set of empirical criteria which separate com- 
posite from "pure" systems. Most of the diagnostic power resides in the W K < 10 
A limit. Other properties (small emission line W FC , low excitation and large Lpm) 
serve to confirm the W^-based diagnostic of compositeness and to judge border line 
cases. We predict that many more starburst /Seyfert 2 composites will be uncovered 
by applying these diagnostics to large samples of near UV spectra. 

(11) The tendency of composites with older starbursts to have smaller H/3, Hell and [OIII] 
emission-line equivalent widths was tentatively interpreted as a sign of the evolution of 
circum-nuclear starbursts in Seyfert 2's. Circum-nuclear starbursts which are powerful 
today may fade away and, barring evolution of the active nucleus, appear as "pure" 
Seyfert 2's. By improving the dating techniques and enlarging the databases we hope 
to trace this evolutionary sequence and thereby put strong constraints on physical 
scenarios for the starburst- AGN connection. 
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